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INTRODUCTION: 


The  Sixl  homeoprotein  is  an  important  regulator  of  embryonic  development  where  it 
controls  proliferation,  survival,  migration,  invasion,  and  expansion  of  progenitor  cell  populations^'"^. 
Sixl  expression  is  normally  restricted  to  embryogenesis,  however  Sixl  is  found  to  be  re-expressed 
many  cancers,  including  breast  cancer,  where  it  is  overexpressed  in  50%  of  primary  cancers,  and  a 
striking  90%  of  metastatic  lesions^.  Using  cell  culture  and  mouse  models,  we  and  others  have 
shown  that  Sixl  contributes  to  tumorigenesis  by  not  only  promoting  proliferation  and  survival*’,  but 
also  by  inducing  an  EMT-like  transformation  and  an  increase  in  the  cancer  stem  cell  population"^’’, 
which  has  been  linked  to  metastatic  and  invasive  cancers*.  We  have  demonstrated  that  Sixl 
upregulates  TGF[3  signaling,  which  is  required  for  the  Sixl-induced  EMT-like  transformation  and 
increase  in  cancer  stem  cell  populations’.  Interestingly,  we  have  also  demonstrated  that  Sixl  is 
capable  of  switching  TGPp  from  a  tumor  suppressor  to  a  tumor  promoter^,  however  the  mechanism 
of  this  switch  is  still  poorly  understood.  Currently,  we  are  investigating  miRNAs  that  may  be 
involved  in  this  Sixl-induced  metastatic  process.  Using  miRNA  microarray  analysis  and  further 
validation  studies  we  have  identified  a  cluster  of  miRNAs  that  are  regulated  by  Sixl:  the  miR-106b- 
25  cluster.  Located  within  the  intron  of  the  MCM7  gene,  this  cluster  contains  three  miRNAs:  miR- 
106b,  miR-93,  and  miR-25.  Previous  research  suggests  that  these  miRNA  have  pro-oncogenic 
properties,  as  well  as  inhibiting  the  growth  suppressive  functions  of  TGPp  signaling Due  to 
these  findings,  and  other  preliminary  data  in  our  lab,  our  aim  in  this  proposal  is  to  examine  the 
Sixl-upregulated  miR106b-25  cluster  and  determine  if  these  miRNA  may  be  the  mechanism  by 
which  Sixl  switches  TGPp  signaling  from  a  tumor  suppressor  to  a  tumor  promoter.  In  addition  we 
also  sought  to  determine  if  these  miRNA  alone  are  capable  of  activating  TGPp  signaling,  and  if  they 
are  necessary  and/or  sufficient  for  Sixl-induced  metastatic  phenotypes.  The  result  of  most  of  this 
work  is  in  our  now  published  manuscript  in  the  journal  Oncogene  in  2012  (attached)'^.  Since  this 
time,  we  have  expanded  on  this  work,  and  have  continued  research  in  the  area  of  miR-106b-25 
induced  metastatic  phenotypes. 


BODY: 

Aiml:  To  determine  if  the  miR106b-25  cluster  contributes  to  the  activation  of  the  TGPp 
pathway  by  Sixl  (month  1-7) 

(A)  Determine  if  the  miR-106b-25  miRNAs  target  endogenous  Smad6  and  Smad7  (month  1-2) 

This  subaim  is  complete,  and  has  been  published  in  our  recent  manuscript*^.  From 
bioinformatic  analysis  we  have  identified  that  all  three  miRNAs  in  the  miR106b-25  cluster  have 
predicted  binding  sites  in  both  the  Smad6  and  Smad7  gene.  These  genes  are  inhibitory  Smads  of 
the  TGPp  pathway,  therefore,  by  repressing  these  genes  the  miRNAs  may  play  a  role  in  activating 
TGPp  signaling*^.  Upon  further  analysis  of  the  function  of  these  genes,  we  decided  to  focus 
primarily  on  Smad7,  since  this  gene  is  intimately  involved  in  the  TGPp  pathway,  while  Smad6  is 
primarily  involved  in  BMP  signaling*"*.  Since  we  know  that  TGPp  signaling  is  critical  for  Sixl- 
induced  metastatic  phenotypes,  we  have  chosen  to  only  look  at  Smad7  for  this  project.  We  have 
first  shown  that  Smad7  mRNA  is  decreased  in  Sixl  overexpressing  cells  (which  is  also  present  at 
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the  protein  level.  When  the  3’UTR  of  Srnad?  was  cloned  into  a  luciferase  reporter,  we  could  also 
identify  repression  of  this  reporter  (targeting  of  the  3’UTR)  in  MCF7  cells  overexpressing  the 
miR106b-25  cluster  (MCF7-Cluster)  versus  an  empty  vector  (MCF7-EV)  and  non- silencing  control 
(MCF7-NS).  Additionally,  we’ve  determined  that  Smad7  protein  is  also  diminished  in  MCF7- 
Cluster  cells  by  western  blot  analysis.  Therefore,  we  have  conclusively  shown  that  the  miR106b-25 
cluster  targets  Smad7  for  repression.  Please  see  Figure  3  in  our  Oncogene  manuscript  for  data 
related  to  this  aim'^. 


B)  Determine  if  introduction  of  miR-106b-25  cluster  miRNAs  into  wild  type  MCF7  cells 
can  induce  TGF-beta  signaling  (month  2-7) 

This  subaim  has  been  completed,  and  is  published  in  our  recent  manuscript*^.  To  address 
this  aim,  we  overexpressed  the  miR106b-25  cluster  both  transiently  and  stably  into  MCF7  cells. 
Upon  overexpression  we  have  demonstrated  that  TGF [3-receptor  type  I  (TpRI)  protein  is 
upregulated  with  miR106b-25  overexpression  (Fig4a-b  in  manuscript*^).  We  have  recently 
published  data  demonstrating  that  TpRI  overexpression  by  Sixl  is  sufficient  for  TGFp  activation^. 
Therefore,  we  propose  that  the  miR106b-25  miRNAs,  at  least  in  part,  contributes  to  the 
upregulation  of  TpRI  protein  induced  by  Sixl  overexpression.  In  addition,  we  also  demonstrate  an 
upregulation  of  phosphorylated-Smad3  (p-Smad3)  in  MCF7-Cluster  cells,  which  is  a  molecular 
marker  for  activated  TGFp  signaling  (Fig4b  in  manuscript*^).  To  strengthen  this  data,  we  have  also 
recently  added  data  showing  that  the  miR106b-25  miRNAs  are  also  sufficient  to  upregulate  many 
genes  that  are  downstream  transcriptional  targets  of  TGFp  signaling,  as  shown  by  quantitative  real¬ 
time  PCR.  Thus,  we  have  also  now  demonstrated  that  these  miRNA  can  activate  TGFp  signaling 
down  to  the  transcriptional  level  (Fig4e  in  manuscript*^). 

Lastly,  we  have  also  conducted  a  microarray  screen  with  our  MCF7-NS  and  MCF7-Cluster 
cells  in  order  to  determine  if  the  miR106b-25  cluster  has  any  impact  on  the  TGFp  response 
signature  (TpRS)  *.  After  analysis,  we  have  shown  that  miR106b-25  overexpression  has  a 
differential  response  signature  than  the  MCF-NS  control  cells  (Supplementary  Figure  7  in 
manuscript*^).  Our  analysis  shown  in  this  figure  is  depicting  only  a  0.5  fold  change,  however,  we 
believe  that  gene  changes  due  to  miRNA  overexpression  are  not  as  strong  as  protein  changes. 
Together,  these  data  demonstrate  that  the  miR106b-25  cluster  can  induce  TGFp  signaling. 

(C)  Determine  if  knock  down  of  miR-106b-25  miRNAs  in  MCF7-Sixl  cells  will  reverse  Six-1 
induced  increases  in  TGF-beta  signaling  (month  2-7) 

This  subaim  has  been  completed,  and  has  been  published  in  our  recent  manuscript*^.  To 
address  this  aim,  we  looked  at  the  same  molecular  markers  as  subaim  B,  which  are  TpRI  protein 
and  p-Smad3.  From  previous  data  in  our  lab,  we  know  that  Sixl  overexpression  will  upregulate 
TpRI  protein  as  well  as  increase  the  expression  p-Smad3,  both  indications  of  increased  TGFp 
signaling.  Therefore,  to  determine  if  the  miRNAs  will  reverse  Sixl-  induced  increases  in  TGFp 
signaling,  we  inhibited  the  miRNAs  (individually,  and  together)  in  the  MCF7-Sixl  cell  line.  For 
miRNA  inhibition,  we  used  the  miRzip  lentiviral  system  from  Open  Biosystems.  We  prepared 
stable  miRzip  cell  lines  in  our  MCF7-Ctrl  and  MCF7-Sixl  cell  lines.  For  a  control,  we  used  a 
miRzip-Scramble  construct,  which  does  not  inhibit  any  miRNA.  Western  blot  of  these  cell  lines 
demonstrates  that  inhibition  of  miR-93  and  the  entire  cluster,  will  reverse  the  TpRI  upregulation 
seen  with  Sixl  overexpressing  cells  (Fig  4C  in  manuscript*^).  Additional  western  blots  with  these 
lysates  indicate  that  knockdown  of  miR-106b,  miR-93,  as  well  as  the  entire  Cluster,  will  reverse  the 
p-Smad3  upregulation  as  seen  with  the  Sixl  overexpressing  MCF7  cells  with  SCR  control  miRzip 
(MCF7-Sixl-SCR)  (Fig  4D  in  manuscript*^).  Therefore,  based  on  these  molecular  markers  for 
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TGPp  signaling,  we  conclude  that  we  do,  indeed,  see  a  reversal  of  Sixl-  induced  increased  in  TGPp 
signaling. 

Aim2:  To  determine  if  the  Sixl-regulated  miR-106b-25  cluster  miRNAs  mediate  the  pro¬ 
metastatic  EMT  functions  of  TGF-beta  while  inhibiting  its  tumor  suppressive 
functions  (month  8-26) 

(A)  Determine  if  addition  of  miR-106b-25  miRNAs  in  wild  type  MCF7  will  recapitulate  the 
EMT  phenotype  of  Sixl  overexpression(month  8-12) 

This  subaim  is  complete  and  has  been  published  in  our  recent  manuscript'^.  The  EMT 
phenotype  as  seen  with  Sixl  overexpression  includes  a  loss  of  E-cadherin  and  P-catenin  from  the 
membrane  (insoluble  fraction),  and  a  gain  of  these  proteins  in  the  soluble  fraction  of  the  celF. 
Additionally,  we  also  report  an  increase  in  P-catenin  transcriptional  activity,  via  a  Topflash 
reporter,  as  well  as  a  decrease  in  adhesion  to  cell  matrix  proteins’.  To  address  whether  or  not  the 
miR106b-25  cluster  can  recapitulate  these  EMT  phenotypes,  we  performed  the  same  assays  with 
MCF7-Cluster  and  MCF7-NS  cells  described  above.  Fractionation  of  the  soluble  and  insoluble 
cell  lysates  shows  a  decrease  of  E-cadherin  and  P-catenin  from  the  membrane  (Insoluble 
fraction),  however  does  not  show  an  increase  of  these  proteins  in  the  soluble  fraction  (Figure  5a 
from  manuscript'^).  Importantly,  however,  we  also  demonstrate  that  miR106b-25  cluster 
overexpression  will  increase  P-catenin  transcriptional  activity  as  measured  by  the  Topflash 
luciferase  reporter  (Fig  5b  from  manuscript'^).  In  addition  to  these  assays,  we  have  also 
performed  immunofluorescence  with  E-cadherin  and  P-catenin,  and  have  demonstrated  a  loss  of 
both  of  these  proteins  from  the  membrane  of  MCF7-Cluster  cells  (Fig5d  from  manuscript'^).  In 
addition  to  this  we  have  also  explored  more  genes  involved  in  EMT  in  these  cells  by  conducting 
an  EMT  gene  PCR  array  (SABiosciences).  Through  this  array,  we  have  now  reported  that  the 
miR106b-25  miRNAs  will  upregulate  many  EMT  related  genes  in  MCF7  cells  (Fig5e).  Lastly, 
we  have  completed  an  adhesion  assay  on  our  MCF7  cells,  now  with  an  additional  control 
(MCF7-EV  for  empty  vector  control).  These  results  now  show  that,  like  Sixl,  miR106b-25 
overexpression  will  result  in  decreased  adhesion  to  Collagen  I,  Collagen  IV,  and  Fibronectin 
(Figure  5f  from  manuscrip'^).  These  phenotypes  suggest  that  the  miRI06b-25  cluster  can  induce 
many  features  of  EMT,  similar  to  Sixl  overexpression  in  MCF7  cells. 

(B)  Determine  if  knock  down  of  miR106b-25  miRNAs  in  MCF7-Sixl  cells 
reverses  the  EMT  induced  by  Sixl  (month  12-16) 

This  subaim  has  been  partially  completed,  and  results  represented  here  have  been 
published  in  our  recent  manuscript'^.  To  determine  if  the  miR106b-25  cluster  is  necessary  for 
Sixl -induced  EMT  phenotypes,  our  plan  was  to  do  the  same  experiments  as  discussed  in  Aim2A, 
but  now  with  an  inhibition  of  the  miRNAs  in  Sixl  overexpressing  cells.  This  task  proved  to  be 
more  difficult  than  miRNA  overexpression  since  the  technology  for  stable  miRNA  inhibition 
was  not  frilly  developed  when  we  started  this  project.  However,  we  were  able  to  start 
experiments  with  transient  oligonucleotide  inhibitors.  Namely,  we  were  able  to  show  that 
inhibition  (via  transient  oligo  inhibitors)  of  all  three  miRNAs  together  in  MCF7-Sixl  cells  will 
decrease  the  Topflash  activity  back  to  MCF7-control  levels  (Figure  5c  from  manuscript'^). 

Because  this  story  has  now  been  published,  we  are  no  longer  working  on  doing  anymore  in  vitro 
assays  with  EMT  features,  we  have  instead  focused  our  efforts  on  our  in  vivo  metastasis 
phenotypes,  which  provide  more  information  on  the  necessity  of  the  miRNAs  than  our  in  vitro 
EMT  assays. 
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(C)  Determine  by  an  in  vivo  model,  if  the  growth  inhibitory  effects  of  TGF-beta 
can  be  restored  in  MCF7-Sixl  cells  by  knock  down  of  miR-106b-25  cluster 
miRNAs  (month  16-26) 

This  subaim  has  been  partially  completed.  Interestingly,  Sixl  overexpressing  cells 
exhibit  a  loss  of  TGPp  growth  inhibition  to  which  the  mechanism  is  unknown.  When  cancer 
cells  switch  TGPp  signaling  from  a  tumor  suppressor  to  a  tumor  promoter,  this  is  known  as  the 
TGPp  paradox.  Previous  research  has  described  the  miR106b-25  cluster  as  having  the  ability  to 
switch  TGPp  from  a  tumor  suppressor  to  a  tumor  promoter  via  inhibition  of  the  cell  cycle 
inhibitor,  p21,  and  the  pro-apoptotic  factor,  Bim*\  Therefore,  we  had  speculated  that  the  Sixl- 
upregulated  miR106b-25  cluster  might  be  a  mechanism  by  which  Sixl  overexpressing  cells 
escape  TGPp  growth  inhibition  (growth  suppressive  functions  of  TGPP).  To  address  this 
question,  we  had  previous  knowledge  that  MCP7  (as  well  as  most  cells  in  vitro)  do  not  respond 
to  TGPp  growth  inhibition  in  vitro.  However,  we  do  know  that  we  are  able  to  see  these  effects  in 
vivo.  Therefore,  we  used  and  in  vivo  system  to  answer  this  question.  Unfortunately,  as 
mentioned  in  Aim2B,  we  originally  had  problems  finding  a  stable  knockdown  of  the  miRNAs,  so 
that  we  may  use  them  in  vivo.  Therefore,  we  chose  a  system  in  which  we  could  stably 
overexpress  the  miRNA,  to  answer  the  same  question.  Previous  data  in  the  lab  has  shown  that 
when  MCP7  cells  are  overexpressed  with  TpRI  protein,  and  then  injected  in  vivo,  that  these  cells 
will  turn  on  TGPp  signaling  and  have  a  growth  suppressive  effect  versus  the  GFP  control. 
Therefore,  we  decided  to  use  our  stable  lentiviral  miRNA  overexpression  constructs 
(miRexpress,  Open  Biosystems)  in  order  to  overexpress  the  miRNA  in  the  context  of  T[3RI 
overexpression.  We  speculated  that  if  the  miRNAs  alone  are  able  to  overcome  TGPp  growth 
inhibition,  that  we  would  see  an  increase  in  the  tumor  growth  in  MCF7-TBRI-Cluster  cells 
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Figure  1:  Cluster  overexpression  in  MCF7-TBRI  cells  does  not  reverse  TGFp  growth  inhibition  in  vivo  as 
measured  by  tumor  growth  (A)  MCF7  cells  overexpressing  either  GFP  control  or  TBRI,  and  overexpressing  NS 
eontrol  or  the  miR106b-25  eluster,  were  injeeted  into  the  mammary  fat  pad  of  NOD-SCID  miee.  Tumor 
measurements  over  time  are  shown  for  eaeh  eell  group  by  ealiper  measurement  (B)  same  tumors  as  in  (A)  but  now 
measured  with  I  VIS  imaging. 

versus  MCF7-TBRI-NS  cells.  In  this  experiment,  we  injected  cells  into  the  4^^  mammary  fat  pad 
of  NOD-SCID  mice,  and  allowed  tumors  to  grow.  We  measured  tumor  size  over  time  using 
caliper  measurements.  Additionally,  we  also  engineered  these  cells  to  have  a  luciferase  reporter 
so  that  we  may  also  measure  tumor  size  by  IVIS  imaging,  after  delivering  luciferase  substrate  to 
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the  mice  via  IP  injection.  We  used  both  methods  to  measure  tumor  volume  over  time,  and 
unfortunately  we  did  not  see  a  reversal  of  tumor  growth  in  the  MCF7-T[3RI-Cluster  cells  versus 
the  MCF7-T[3RI-NS  cells  as  measured  by  a  caliper  (Fig  lA).  The  IVIS  measurements,  which  are 
measuring  the  luciferase  signal  of  the  tumor  cells,  were  inconclusive,  as  these  measurements  did 
not  correlate  with  the  actual  tumor  size  (Fig  IB).  We  believe  that  the  tumors  in  these  mice  were 

too  large  to  measure  by  luciferase  signal  because 
they  contained  necrotic  regions.  This  hypothesis 
was  verified  with  tumor  dissections  to  identify 
many  necrotic  regions  in  several  of  the  tumors. 
From  these  experiments,  we  did  not  identify  the 
miR106b-25  cluster  as  being  sufficient  for  TGF[3 
growth  inhibition  in  vivo.  However,  we  believe 
that  the  upregulation  of  this  cluster  by  Sixl  may 
still  be  necessary  for  Sixl’s  ability  to  evade  TGF[3 
growth  inhibition. 

Since  our  last  report,  we  have  additional 
evidence  of  the  role  of  the  miR106b-25  miRNAs  in 
TGFp  growth  suppression,  and  we  have  shown  this 
by  utilizing  cell  lines  that  do  indeed  respond  to 
TGFp  in  a  growth  suppressive  manner,  in  culture. 
One  of  these  cell  lines  is  the  normal  murine 
mammary  gland  cell  line  (NMuMGs).  In  this 
system  we  have  shown  that  when  we  overexpress 
Sixl  in  these  cells  (NMuMG-Sixl),  we  have  a 
decreased  sensitivity  to  TGF-P  growth  inhibition  as 
compared  to  control  (NMuMG-Ctrl),  as  seen  by  a 
3H-Thymidine  incorporation  assay  with  increasing 
TGF-P  concentrations  (Figure  2,  top).  In  the 
NMuMG-Sixl  cells,  we  have  now  utilized  locked 
nucleic  acid  (LNA)  inhibitors  (provided  by 
miRagen®)  to  knock  down  the  miR106b-25 
miRNAs  in  the  context  of  Sixl .  When  we  knock 
down  all  three  miRNA  (miR-106b,93,  and  25)  in 
the  NMuMG-Sixl  cells,  we  observe  a  decreased 
sensitivity  to  TGF-P  growth  inhibition  (Figure  2, 
bottom),  suggesting  that  these  miRNA  do  indeed 
play  a  role  in  Sixl-mediated  TGF-P  growth 
suppression. 

To  strengthen  this  argument  we  have  also 
identified  a  breast  cancer  cell  line  which  contains 
both  high  levels  of  Sixl  and  the  miR-106b-25 
miRNAs,  the  BT549  cells.  In  a  panel  of  breast 
cancer  cell  lines,  these  cell  lines  exhibit  a  higher 
expression  of  Sixl/miR-106b-25  than  other  cell 
lines  (Figure  3A).  In  this  context,  we  have  now 
also  shown  that  knockdown  of  the  miR-106b-25 
miRNAs  with  LNA  inhibitors  also  decreases  the  sensitivity  of  these  cells  to  TGF-P  growth 
inhibition  (Figure  3B).  Together,  these  additional  in  vitro  results  strengthen  the  hypothesis  that 


NMuMG-GFP  vs.  Sixl:  TGFp  growth  inhibition 


0  0.05  0.5  5  10 

ng/mLTGFpi 


ng/mLTGFpi 

Figure  2.  The  miR-106b-25  miRNAs  are 
required  for  Sixl  to  overcome  TGF-p- 
mediated  growth  suppression,  (a)  Normal 
Murine  Mammary  Gland  cells  (NMuMG) 
overexpressing  Sixl  (NMuMG-Sixl)  or 
expressing  a  control  plasmid  (NMuMG-GFP) 
were  subjected  to  increasing  concentrations  of 
TGF-pl  (added  to  the  media).  After  48  hours 
cells  were  assessed  for  their  proliferative 
potential  through  3H-Thymidine  incorporation. 
Each  cell  line  was  plotted  as  %  growth  at 
individual  concentrations  of  TGF-pl  as  compared 
to  vehicle  treated  cells  with  no  TGF-pl  addition. 
NMuMG-GFP  cells  (red)  exhibit  a  decrease  in 
growth  with  TGF-pl  addition,  while  NMuMG- 
Sixl  (blue)  cells  exhibit  resistance  to  growth  as 
compared  to  GFP  control  cells  (b)  miR-106b-25 
cluster  miRNAs  were  inhibited  in  NMuMG-Sixl 
cells  with  LNA-miRNA  inhibitors.  Inhibition  of 
all  three  miRNAs  +  their  family  members 
(106b+25  LNA)  in  NMuMG-Sixl  cells  sensitized 
these  cells  to  TGF-pl  growth  inhibition,  as  read 
bv  MTS  assay,  and  compared  to  NMuMG-Sixl 
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the  miR-106b-25  miRNAs  are  necessary  for  Six  1 -mediated  TGF-[3  growth  suppression,  and 
move  us  closer  to  understanding  the  mechanism  of  the  TGF-[3  switch  in  the  context  of  Sixl 
overexpression. 


a 

50.,  miR-106b  ,5  rniR-93  ,(w  miR-25 


BT549  cells  with  miRNA  LNA  inhibitors 


Figure  3.  miR-106b-25  miRNA  inhibition  sensitizes 
BT549  cells  to  TGF-p  growth  inhibition,  (a)  All  three 
miRNAs  in  the  miR-106b-25  eluster  have  high  endogenous 
expression  in  the  BT549  eell  line,  as  eompared  to  other 
breast  eaneer  eells  (above),  whieh  eorresponds  to  a  high 
endogenous  expression  of  Sixl  in  this  partieular  eell  line 
(bottom).  qRT-PCR  was  performed  for  expression  analysis, 
with  normalization  to  U6  RNA  for  miRNA  expression,  and 
to  PPIB  for  Sixl  expression  (b)  BT549  eells  with  inhibition 
of  miR-106b/93  (106b  Faml),  miR-25  (25),  or  all  together 
(106b  Fam+25),  exhibit  deereased  growth  with  inereasing 
eoneentrations  of  TGF-pi,  as  measured  by  MTS  assay,  and 
eompared  to  BT549  eells  transfeeted  with  a  eontrol  LNA 
(Trune_Ctrl).  Graph  represents  %  growth,  as  measured  in 
Figure  B.l. 


Aim3:  To  determine  the  in  vivo  role  of  the  miR106b-25  cluster  miRNAs  in  Sixl-lnduced 
metastasis  (month  26-36) 

(A)  Determine  if  mlR106b-25  miRNA  upregulatlon  is  sufficient  to  induce  metastasis  in  an 
orthotopic  xenograft  breast  cancer  mouse  model  (month  26-31) 

This  subaim  is  mostly  completed.  To  answer  this  question  our  first  approach  was  to  use  an 
experimental  metastasis  model,  in  which  we  inject  luciferase-taggeed  cells  into  the  left  ventricle 
of  the  heart  of  nude  mice,  and  follow  metastasis  of  cells  through  I  VIS  imaging.  We  have 
previously  published  this  metastasis 
model  and  the  ability  of  MCF7-Sixl  cells 
to  increase  metastasis  over  MCF7-Ctrl 
cells  in  this  model’.  Therefore,  we 
decided  to  see  if  the  MCF7-Cluster  cells 
would  increase  metastasis  over  MCF7- 
NS  cells  in  this  model  system,  alongside 
MCF7-Ctrl  and  MCF7-Sixl  cells  to  be 
able  to  compare  the  metastatic  results. 

To  achieve  this,  we  luciferase-tagged  our 
MCF7-NS  and  MCF7-Cluster  cells,  and 
used  the  MCF7-Ctrl  and  MCF7-Sixl 
luciferase  tagged  cells  from  the  published 
experiment.  We  used  all  4  cell  lines  and 
injected  these  cells  intracardiac  into  nude 


MCF7-CTRL 

MCF7-Six1 

MCF7-NS 

MCF7-Cluster 


Figure  4:  An  experimental  metastasis  model  does  not  show 
increased  metastasis  with  Cluster  overexpression. 

Luciferase  tagged  MCF7  eells  overexpressing  Ctrl,  Sixl,  NS, 
and  Cluster  eonstruets  were  injeeted  intraeardiae  into  nude 
miee.  Metastasis  of  these  eells  were  measured  by  IVIS 
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mice  (15  mice  per  group).  We  were  able  to  measure  metastasis  via  IVIS  imaging  over  time  (Fig 
4).  As  before,  we  were  able  to  see  an  increase  in  metastasis  in  the  MCF7-Sixl  cells  over  time, 
however,  we  did  not  see  any  increase  in  metastasis  with  MCF7-Cluster  cells  over  their  control 
cells  (MCF7-NS)  (Fig  4).  The  experimental  metastasis  model,  however,  is  an  assay  designed  to 
measure  the  late  stages  of  metastasis.  When  cells  are  injected  directly  into  the  blood  stream,  the 
initial  metastatic  stages  of  intravasation,  motility,  and  extravasation,  are  bypassed.  Therefore, 
we  cannot  disclude  the  possibility  that  the  miR-106b-25  cluster  may  be  important  for  the  earlier 
stages  of  metastasis. 

We  have,  however,  conducted  orthotopic  experiments  for  the  purpose  of  discerning 
whether  or  not  the  miR-106b-25  cluster  is  sufficient  to  initiate  tumors  at  limiting  dilutions,  which 
is  a  common  assay  used  for  cancer  stem  cell  analysis.  This  assay  did  not  include  a  metastasis 
element,  but  it  did  demonstrate  that  these  miRNAs  have  an  increased  population  of  cancer  stem 
cells,  similar  to  Sixl,  as  the  mice  that  were  injected  with  MCF7-Cluster  cells  had  a  significant 
increase  in  tumor  formation  over  the  NS  control  cells  injected  (Figure  6C  in  manuscript'^).  This 
assay  not  only  demonstrates  that  the  miR106b-25  cluster  also  increases  the  cancer  stem  cell 
population,  but  that  we  can  demonstrate  an  increase  in  tumor  formation  in  an  orthotopic  model. 

(B)  Determine  if  miR-106b-25  miRNAs  are  necessary  for  Sixl-induced  metastasis  in  an 
orthotopic  xenograft  breast  cancer  mouse  model  (month  31-36) 

This  subaim  remains  unanswered  due  to  unforeseen  complications.  Since  our  last  report, 
we  have  conduced  an  in  vivo  experiment  to  test  the  necessity  of  the  miR-106b-25  miRNAs  in 
Sixl -mediated  metastasis.  Our  approach,  however,  was  not  an  orthotopic  xenograft  model.  We 

chose  to  instead  use  our  experimental  metastasis  model 
described  above.  Our  reasons  for  this  included  the  fact  that 
this  is  a  well  established  model  system  in  our  laboratory  for 
Sixl -mediated  metastasis.  In  order  to  survey  the  necessity 
of  the  miR-106b-25  miRNAs  in  this  system,  we  decided  to 
take  advantage  of  a  collaboration  with  miRagen 
therapeutics,  who  has  provided  us  with  in  vivo  deliverable 
miRNA  inhibitors.  Thus,  miRagen  designed  and 
manufactured  specific  LNA  oligos  that  would  inhibit  the 
miR-106b-25  cluster  miRNAs  for  our  use  in  vivo.  In  our 
experimental  design,  we  again  injected  MCF7-Ctrl  and 
MCF7-Sixl  (which  contain  high  miR-106b-25  miRNAs) 
cells  in  the  left  ventricle  of  nude  mice  at  DayO.  Starting 
Day3,  we  treated  our  experimental  group  of  animals  with 
our  LNA  inhibitors  through  subcutaneous  injections  in  the 
rear  flank  of  the  animals.  Unfortunately,  soon  after  LNA 
delivery,  the  animals  developed  a  reaction  to  the  LNA 
injections.  Bruising  and  dark  lesions  occurred  at  the  site  of 
delivery  (Figure  5).  These  unforeseen  injuries  lead  to  a 
disruption  of  LNA  delivery,  as  our  animal  facility  halted 
drug  treatment  in  the  second  week  of  our  experiment  (LNA 
delivery  originally  planned  bi-weekly).  We  continued  the 
experiment  in  week3,  and  continued  to  treat  the  animals 
with  LNA  inhibitors  and  used  IVIS  imaging  for  monitoring 
metastasis.  Our  final  results  in  metastasis  did  not  reveal  any 
significant  changes  with  Sixl -mediated  metastasis  in  the 


Figure  5.  Unexpected  lesions  in 
mice  as  caused  by  miR-106b  family 
LNA  inhibitor  subcutaneous 
injections.  (Top)  The  picture  reveals 
a  nude  mouse  with  a  lesion  on  its  left 
hind  flank  region  at  a  24-hour  time 
point  following  LNA  subcutaneous 
injection.  (Bottom)  Lesions  continued 
to  heal  over  time,  and  some  lesions 
developed  into  tough  scabbing  around 
1  week  post  LNA  injection. 
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miR  inhibited  group  (MCF7-Sixl-106b+25  Fam)  as  compared  to  control  Sixl  overexpressing 
group  (MCF7-Sixl-SCR)  (Figure  6).  Although,  we  believe  that  disruption  of  LNA  delivery  in 
the  second  week  of  treatment  could  have  been  crucial  in  miR-induced  colonization  of  metastasis. 
Therefore,  we  consider  this  aim  incomplete,  and  plan  to  conduct  this  experiment  differently, 
with  genetic  inhibition  of  miRs  I  cancer  cells  instead  of  systemic  LNA  miR  inhibition. 

Our  results  in  the  LNA  experiment,  however,  bring  about  important  questions  concerning 
the  possible  side-effects  of  miR106b-25  inhibition.  After  pathological  evaluation  of  our  LNA 
specific  lesions,  we  determined  that  the  lesions  were  only  caused  by  specific  inhibition  of  the 
miR-106b/93  family  of  miRNAs.  We  also  have  determined  that  this  was  a  skin  specific  lesion. 

Therefore,  we  believe  that  we  have 
uncovered  possible  crucial 
information  concerning  miR-106b/93, 
and  it’s  role  in  the  skin.  We  believe 
this  is  an  important  avenue  to  follow¬ 
up  with  for  future  considerations  of 
inhibiting  these  miRNA 
therapeutically. 

Importantly,  we  do,  however, 
have  additional  data  suggesting  that 
the  miR-106b-25  miRNA  are  indeed 
necessary  for  Sixl  metastatic 
phenotypes.  We  have  yet  to 
determine  the  actual  impact  on 
metastasis,  however,  we  have 
discovered  that  the  miR-106b-25 
miRNAs  are  necessary  for  Sixl- 
induced  tumor  initiation  (cancer  stem 
cell  phenotype).  We  have 
demonstrated  this  in  our  latest 
manuscript  by  showing  that  MCF7- 
Sixl  cells  with  miR-106b-25 
knockdown  (with  stable  miRZip 
lentiviral  constructs)  have  a  decreased 
ability  to  form  tumors  at  limiting 
dilutions,  as  compared  to  control  cells 
overexpressing  Sixl  (Figure  6D  in  manuscript'^).  Since  the  increase  of  cancer  stem  cells  is  an 
important  property  of  metastasis,  we  consider  this  to  be  strong  evidence  that  the  miR-106b-25 
miRNA,  at  least  in  part,  contribute  to  Sixl -induced  metastasis. 

Key  Research  Accomplishments:  All  accomplishments  from  start  of  funding  are  listed, 
accomplishments  since  last  report  are  in  bold  type. 

•  Published  many  of  these  findings  in  a  manuscript  entitled:  “The  miR-106b-25  cluster 
targets  Smad7,  activates,  TGF-beta  signaling,  and  induces  EMT  and  tumor  initiating  cell 
characteristics  downstream  of  Sixl  in  human  breast  cancer’’  in  the  journal  Oncogene  in 
2012. 

•  Determined  that  the  miR106b-25  cluster  targets  the  Smad7  gene  for  repression 

•  Determined  that  overexpression  of  the  miR106b-25  cluster  is  sufficient  to  upregulate 


Metastasis  as  measured  by  total  flux 


Figure  6.  miR-106b-25  inhibition  through  LNA  delivery 
reveals  no  difference  in  metastasis  in  MCF7-Sixl  injected 
mice.  In  an  experimental  metastasis  model,  miee  injeeted  with 
MCF7-Ctrl  or  MCF7-Sixl  eells  were  treated  with  LNA  miR 
inhibitors  that  were  either  a  serambled  etrl  (SCR),  a 
eombination  of  miR- 106b  family  and  miR-25  family  inhibitors 
(106b+25Fam),  or  saline  only.  The  graph  represents  total  flux 
lueiferase  measurement  in  eaeh  animal  group  over  time.  Total 
flux  values  +/-  standard  deviation  for  eaeh  treatment  group  are 
shown  at  eaeh  time  point.  Data  show  no  differenee  in 
metastatie  burden  of  MCF7-Sixl  miee  treated  with  SCR, 
106b+25Fam,  or  saline  only. 
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XpRI  protein,  p-Smad3,  and  many  TGPp  transcriptional  targets.  We  also  demonstrated  a 
differential  response  in  the  TGPp  response  signature,  indicating  activation  of  the  TGPp 
pathway  by  the  miRNAs  alone 

•  Identified  that  inhibition  of  miRNAs  in  Six  1 -overexpressing  MCP7  cells  will  reverse  the 
activation  of  TGPp  signaling  as  measured  by  decreased  XpRI  protein  and  p-Smad3 

•  Determined  that  overexpression  of  the  miR106b-25  cluster  alone  will  induce  a  loss  of  E- 
Cadherin  and  P-catenin  from  the  membrane,  as  well  as  an  increase  in  P-catenin 
transcriptional  activity.  We  have  also  shown  by  immunofluorescence  a  loss  of  these 
proteins  from  the  membrane  of  MCP7-Cluster  cells.  We  have  further  shown  that  the 
miR106b-25  miRNAs  induce  EMX  through  an  induction  of  many  EMX  related  genes  and 
a  reduction  of  adhesion  to  cell  matrix  proteins.  Xhus,  conclusively  indicating  that  these 
miRNAs  are  sufficient  to  induce  features  of  EMX. 

•  Determined  that  knock  down  of  the  miR106b-25  cluster  is  sufficient  to  reverse  P-catenin 
transcriptional  activity  in  Sixl  overexpressing  MCP7  cells 

•  Identified  that  the  miR106b-25  cluster  alone  is  not  sufficient  to  increase  metastasis  in  an 
experimental  metastasis  model  (late  stage  metastasis) 

•  Determined  that  the  miRI06b-25  miRNAs  are  sufficient  AND  necessary  for  tumor 
initiation  (or  cancer  stem  cell  capacity)  in  an  in  vivo  tumor  initiating  assay. 

•  Now  have  further  evidence  in  support  of  the  necessity  of  the  miR-106b-25  miRNAs 
in  Sixl-mediated  TGF-P  growth  inhibition 

•  Have  determined  that  miR-106b/93  inhibition  leads  to  adverse  side-effects  in  the 
skin  of  nude  mice,  which  is  a  possible  new  avenue  for  a  role  for  these  miRNA  in  the 
skin. 


Reportable  outcomes  for  2013: 

Pirst  author  review  based  on  XGP-P  biology: 

Smith  AL,  Robin  XP,  &  Pord  HE.  Molecular  Pathways:  Xargeting  the  XGP-P  pathway  for 
Cancer  Xherapy.  Clin  Cancer  Res.  18,4514-4521  (2012) 

Research  presented  in  the  form  of  a  poster: 

Smith  AL,  Iwanaga  R,  Drasin  DJ,  Micalizzi  DS,  Vartuli  RE,  Pord  HE.  Xhe  miR-106b-25 
Cluster  as  a  Potential  Xarget  for  Sixl-mediated  Breast  Cancer  Metastasis,  and  a  Molecular 
Marker  of  the  XGPP  Paradox.  EORXC-NCI-AACR  Molecular  Xargets  and  Cancer  Xherapeutics, 
Dublin,  Ireland,  November  6-9,  2012  See  appendix  for  poster  abstract. 

•  Bolie  graduate  scholarship  fund  travel  award  (from  Molecular  Biology  Program  at 
UCDenver  AMC)  awarded  for  travel  to  this  conference 

Research  presented  in  the  form  of  a  thesis  defense: 

•  March  E*,  2013:  “Xhe  role  of  the  miR-106b-25  miRNAs  in  Sixl-mediated  breast  cancer 
progression  and  metastasis” 

•  phD  degree  in  Molecular  Biology  obtained 

•  post-doctoral  positions  offered  based  research  in  this  project  (no  decisions  made  yet) 
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Conclusions: 


The  results  from  studies  up  to  this  point  have  shown  that  the  miRNAs  (miR-106b,  miR-93,  and 
miR-25)  all  have  an  impact  on  Sixl  induced  tumorigenesis.  We  have  demonstrated  that  the 
miR106b-25  cluster  can  repress  the  TGF[3  inhibitor  Smad7,  thus  having  an  effect  on  activation  of 
the  TGPp  pathway.  Indeed,  we  have  demonstrated  that  the  miR106b-25  cluster  is  both  sufficient 
and  necessary  for  activation  of  TGPp  signaling.  We  have  conclusively  shown  that  these  miRNA 
are  sufficient  for  EMT  induction,  similar  to  Sixl,  through  a  loss  of  both  E-cadherin  and  P-catenin 
from  the  membrane,  a  loss  of  cell  matrix  adhesion,  and  an  upregulation  of  several  genes  involved 
in  EMT.  This  work  has  also  lead  to  our  increased  understanding  of  a  role  for  the  miR-106b-25 
miRNAs  in  cancer  stem  cell  biology.  We  have  now  demonstrated  that  the  miR-106b-25 
miRNAs  are  both  necessary  and  sufficient  for  Sixl -induced  cancer  stem  cell  phenotypes.  We  do 
no  not  yet  know  if  these  phenotypes  are  through  the  miR-induced  increase  in  TGF-P  signaling. 
Importantly,  we  have  made  continued  progress  this  year  in  the  attempts  to  determine  the  role  of 
the  miRs  in  Sixl-induced  TGPp  growth  suppression.  We  now  have  further  evidence  in  two 
additional  cell  lines  confirming  that  the  miR-106b-25  miRNAs  are  important  for  the  TGF-P 
switch  downstream  of  Sixl. 

We  were  not  able  to  demonstrate  that  the  miR106b-25  cluster  is  sufficient  for  increased 
metastasis  through  an  experimental  metastasis  model.  However,  it  is  possible  that  these 
miRNAs  may  only  be  necessary  and  not  sufficient  for  metastasis.  When  attempting  to  determine 
the  necessity  of  the  miRs  in  Sixl-induced  metastasis,  we  uncovered  a  possible  side-effect  of 
miRNA  inhibition  using  ENA  miRNA  inhibitors  delivered  sub-cutaneously.  We  believe  that  this 
disturbance  caused  our  results  in  this  experiment  to  be  interpretable,  and  therefore  we  plan  to 
conduct  this  experiment  again,  but  with  genetic  downdown  of  the  miRNAs  instead  of  systemic 
inhibition  with  LNAs.  Since  we  have  shown  that  these  miRNA  are  important  for  TGPp 
signaling,  features  of  EMT,  as  well  as  cancer  stem  cell  phenotypes,  we  anticipate  that  we  will  see 
an  effect  on  Sixl-induced  metastasis  when  inhibiting  the  miRNAs  through  a  different 
mechanism. 
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Appendices: 

Poster  abstract  for  EORTC-NCI-AACR  Molecular  Targets  and  Cancer  Therapeutics,  Dublin, 
Ireland,  November  6-9,  2012: 

The  miR-106b-25  Cluster  as  a  potential  target  for  Sixl -mediated  breast  cancer  metastasis,  and  a 
molecular  marker  of  the  TGFp  paradox. 

Previous  work  in  our  laboratory  has  demonstrated  a  critical  role  for  the  homeoprotein 
Sixl  in  breast  cancer  metastasis.  Our  data  has  shown  that  Sixl  controls  metastatic  spread 
through  an  upregulation  of  TGFp  signaling,  leading  to  an  EMT-like  transformation  and  an 
increase  in  the  tumor  initiating  cell  (TIC)  population.  TGFp  is  a  multifunctional  cytokine,  which 
can  function  as  either  a  tumor  suppressor  or  a  tumor  promoter,  making  this  pathway  a  difficult 
target  in  cancer  therapy.  Additionally,  our  data  shows  that  the  presence  of  Sixl  in  cells  will 
implement  the  switch  of  TGFp  from  tumor  suppressive  to  tumor  promotional. 

To  further  understand  the  mechanism  by  which  Sixl  mediates  metastasis  through  TGFp 
signaling,  we  explored  miRNAs  that  may  be  involved  in  this  process.  Through  a  miRNA 
microarray  and  further  validation  studies,  we  found  that  Sixl  upregulates  a  cluster  of  three 
miRNAs,  the  miR-106b-25  cluster,  which  is  comprised  of  miR-106b,  miR-93,  and  miR25. 
Interestingly,  this  cluster  has  already  been  implicated  in  the  impairment  of  TGFP  growth 
suppression  through  repression  of  the  cell  cycle  inhibitor,  p21,  and  pro-apoptotic  factor,  Bim. 
Therefore  suggesting  that  upregulation  of  these  miRNA  by  Sixl  may  help  mediate  the  switch  in 
TGFp  signaling.  Indeed,  preliminary  results  in  NMuMG  (normal  mouse  mammary  gland)  cells 
and  BT549  breast  cancer  cells  demonstrates  that  inhibition  of  the  miR-106b-25  miRNAs  in  the 
context  of  Sixl  overexpression  restores  TGFP-mediated  growth  inhibition.  In  addition  to  its  role 
in  TGFp  growth  suppression,  we  have  also  recently  identified  that  this  same  cluster  of  miRNA  is 
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capable  of  activating  the  tumor  promotional  functions  of  TGPp.  We  have  shown  that  these 
miRNA  bind  and  repress  the  TGPp  inhibitor,  Smad7,  and  activate  downstream  TGPp  signaling 
through  upregulation  of  p-Smad3  and  TGPp  transcriptional  targets.  Upon  further  exploration, 
we  also  find  that  the  miR-106b-25  cluster  is  sufficient  to  induce  EMT-like  features  as  well  as  an 
increase  in  the  TIC  population,  and  that  these  miRNA  are  required  downstream  of  Sixl  to  induce 
these  phenotypes.  Lastly,  we  also  demonstrate  a  significant  correlation  of  Sixl  and  miR-106b 
expression  in  human  breast  cancers,  and  further  show  that  high  expression  of  miR-106b  or  miR- 
93  correlates  to  shortened  time  to  relapse  in  breast  cancer  patients.  All  together,  we  find  that  the 
Sixl -mediated  increase  in  the  miR-106b-25  cluster  may  not  only  be  involved  in  facilitating  the 
oncogenic  switch  of  TGPp,  but  also  aids  in  the  activation  of  this  pathway,  thus  contributing  to 
the  pro-metastatic  phenotypes  found  with  Sixl  overexpression.  Because  it  is  difficult  to  target 
transcription  factors  like  Sixl  directly,  our  data  points  to  a  new  therapeutic  target  for  Sixl 
expressing  breast  cancers  (miR-106b-25),  which  may  also  be  used  to  predict  which  patients 
should  benefit,  rather  than  be  harmed,  by  TGPp  inhibitors  which  are  currently  in  clinical  trials. 
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ORIGINAL  ARTICLE 

The  miR-106b-25  cluster  targets  Smad7,  activates  TGF-P 
signaling,  and  induces  EMT  and  tumor  initiating  cell 
characteristics  downstream  of  Sixl  in  human  breast  cancer 

AL  Smith\  R  Iwanaga^  DJ  Drasin’,  DS  Micalizzi’'^  RLVartuli’,  A-C  Tan"  and  HL 


The  role  of  TGF-p  signaling  in  tumorigenesis  is  paradoxical:  it  can  be  tumor  suppressive  or  tumor  promotional,  depending  on 
context.  The  metastatic  regulator,  Sixl,  was  recently  shown  to  mediate  this  switch,  providing  a  novel  means  to  explain  this 
elusive  TGF-p  paradox'.  Flerein,  we  identify  a  mechanism  by  which  Sixl  activates  the  tumor  promotional  arm  of  TGF-p 
signaling,  via  its  ability  to  upregulate  the  miR-106b-25  microRNA  cluster,  and  further  identify  a  novel  function  for  this  cluster 
of  microRNAs.  Although  expression  of  the  miR-106b-25  cluster  is  known  to  overcome  TGF-p-mediated  growth  suppression  via 
targeting  p21  and  BIM,  we  demonstrate  for  the  first  time  that  this  same  cluster  can  additionally  target  the  inhibitory  Smad7 
protein,  resulting  in  increased  levels  of  the  TGF-p  type  I  receptor  and  downstream  activation  of  TGF-p  signaling.  We  further 
show  that  the  miR-106b-25  cluster  is  sufficient  to  induce  an  epithelial-to-mesenchymal  transition  and  a  tumor  initiating  cell 
phenotype,  and  that  it  is  required  downstream  of  Sixl  to  induce  these  phenotypes.  Finally,  we  demonstrate  a  significant 
correlation  between  miR-106b,  Sixl,  and  activated  TGF-p  signaling  in  human  breast  cancers,  and  further  show  that  high  levels 
of  miR-106b  and  miR-93  in  breast  tumors  significantly  predicts  shortened  time  to  relapse.  These  findings  expand  the  spectrum 
of  oncogenic  functions  of  miR-106b-25,  and  may  provide  a  novel  molecular  explanation,  through  the  Sixl  regulated  miR-106b-25 
cluster,  by  which  TGF-p  signaling  shifts  from  tumor  suppressive  to  tumor  promoting. 

Oncogene  advance  online  publication,  30  January  2012;  doi:10.1038/onc.2012.1 1 
Keywords:  Sixl;  miRNA;  TGF-P;  epithelial-to-mesenchymal  transition;  tumor-initiating  cell 


INTRODUCTION 

TGF-p  signaling  has  a  critical  and  dual  role  in  breast  tumorigenesis. 
In  normal  epithelium  as  well  as  in  early  tumorigenic  lesions,  TGF-p 
has  a  tumor  suppressive  role  through  its  ability  to  induce  growth 
inhibition.  As  cancer  progresses,  however,  tumor  cells  become 
resistant  to  TGF-p-mediated  growth  inhibition,  and  instead  TGF-p 
promotes  tumor  progression  and  metastasis,  likely  in  part  through 
its  promotion  of  an  epithelial-to-mesenchymal  transition  (EMT).^'^ 
The  ability  of  TGF-p  to  switch  cells  from  tumor  suppressive  to 
tumor  promotional  has  been  coined  the  TGF-p  paradox.^  The 
mechanism  of  this  phenomenon  is  not  well  understood,  but 
remains  an  important  area  of  research. 

Recently,  our  laboratory  identified  the  developmental  homeotic 
transcription  factor,  Sixl,  as  an  important  mediator  of  breast 
cancer  progression  and  metastasis.^“^  We  have  found  that  Sixl  is 
dependent  on  upregulation  of  TGF-p  signaling  to  induce  an  EMT,^ 
an  increase  in  tumor-initiating  cell  (TIC)  characteristics^  (and  data 
not  shown),  and  to  induce  late  stage  metastasis.^  Interestingly, 
using  an  experimental  metastasis  model,  we  demonstrated  that 
Sixl  not  only  upregulates  TGF-p  signaling,  but  that  it  also  switches 
the  pathway  from  tumor  suppressive  to  tumor  promotional.^'^ 
Flowever,  until  now,  the  mechanism  by  which  Sixl  accomplishes 
this  switch  was  unknown. 

Flerein,  we  have  identified  a  cluster  of  miRNAs,  miR-106b-25, 
that  is  upregulated  by  the  homeoprotein  Sixl.  This  cluster 


is  highly  conserved  in  vertebrates,  and  consists  of  three  miRNAs, 
miR-106b,  miR-93  and  miR-25,  which  all  reside  in  the  13th  intron 
of  the  MCMy  gene  (Chr7).  Two  paralogs  of  this  cluster  exist,  the 
miR-17-92  cluster  (Chr13,  c130RF25  gene),  and  the  miR-106a-363 
cluster  (ChrX).  The  miR-17-92  cluster  has  received  considerable 
attention  as  a  pro-oncogenic  cluster,  and  was  the  first  cluster  of 
miRNA  identified  to  cooperatively  act  as  an  oncogene.^®  Interest¬ 
ingly,  the  miR-106b-25  cluster  of  miRNA  is  overexpressed  in 
several  cancers,^ ^  and  like  the  miR-17-92  paralog,  is  pro- 
oncogenic.^^  Previous  reports  have  demonstrated  the  ability  of 
these  miRNA  to  inhibit  the  growth  suppressive  functions  of  TGF-p 
signaling  through  repression  of  downstream  mediators  p21  and 
Bim.^^  In  this  study,  we  show  for  the  first  time  that  this  cluster  of 
miRNA  can  target  the  TGF-p  inhibitor  Smad-7  and  also  activate  the 
TGF-p-signaling  pathway,  providing  a  novel  mechanism  by  which 
Sixl  overexpression  can  mediate  the  switch  in  TGF-p  signaling 
from  tumor  suppressive  to  tumor  promotional.  We  also  demon¬ 
strate  that  overexpression  of  the  miR-106b-25  miRNAs  is  sufficient 
to  induce  characteristics  of  EMT  and  TICs,  and  that  this  cluster  is 
necessary  for  the  ability  of  Sixl  to  mediate  these  phenotypes. 
Finally,  we  demonstrate  a  significant  correlation  between  miR- 
106b  expression  and  both  Sixl  and  activated  TGF-p  signaling  in 
human  breast  cancer  tissues,  and  further  show  that  high 
expression  levels  of  miR-106b  and  miR-93  together  in  early- 
invasive  breast  tumors  can  significantly  predict  a  shorter  time  to 
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relapse  for  these  patients.  Taken  together,  our  results  significantly 
expand  the  tumorigenic  functions  of  the  miR-106b-25  cluster. 
In  addition,  we  demonstrate  a  critical  role  for  this  cluster  in 
mediating  not  only  the  pro-tumorigenic  functions  of  Sixl,  but 
also  provide  a  possible  mechanism  by  which  Sixl  overcomes 
TGF-p-mediated  growth  suppression,  whereas  simultaneously 
activating  the  pro-metastatic  arm  of  the  pathway. 


RESULTS 

Sixl  regulates  the  miR-106b-25  Cluster  of  miRNAs 
Previous  studies  have  demonstrated  substantial  cross-talk 
between  miRNAs  and  homeobox  genes.^^'^^  We  therefore  asked 
whether  the  Sixl  homeoprotein  might  regulate  miRNAs  to 
mediate  its  tumorigenic  and  metastatic  phenotypes.  miRNA 
microarray  analysis  on  RNA  isolated  from  IVICF?  breast  cancer 
cells  overexpressing  Sixl  (MCF7-Six1)  and  control  cells  (MCF7-Ctrl) 
led  to  the  identification  of  several  miRNAs  that  were  differentially 
expressed  in  a  statistically  significant  manner  between  the  two 
groups  (Figure  la).  Interestingly,  we  identified  two  miRNAs,  miR- 
106b  and  miR-25,  that  were  upregulated  in  response  to  Sixl 
overexpression  (Figure  1  a),  and  that  belong  to  a  cluster  of  miRNAs, 
which  also  includes  miR-93,  and  reside  in  the  13th  intron  of  the 
MCM7  gene  (Figure  1b).  These  miRNA  have  previously  been 
implicated  as  a  pro-oncogenic  cluster  of  miRNAs.^^'^^'^^  To  validate 
our  microarray  results,  we  performed  quantitative  real-time 
reverse  transcriptase  PCR  (qRT-PCR)  on  an  independent  set  of 
RNA  isolated  from  MCF7-Ctrl  and  MCF7-Six1  cells,  demonstrating 
that  all  three  miRNA  within  the  cluster  are  overexpressed  2-3  fold 
in  MCF7-Six1  cells  as  compared  with  MCF7-Ctrl  cells  (Figure  2a).  In 
addition,  siRNA  knockdown  of  Sixl  in  21 PT  cells  (Supplementary 
Figure  1),  which  contain  high  levels  of  Sixl  endogenously,^ 
resulted  in  a  clear  decrease  in  all  three  miRNAs,  confirming  that 
endogenous  Sixl  regulates  the  miR-106b-25  cluster  (Figure  2b). 
Finally,  to  examine  whether  Sixl  could  regulate  the  miR-106b-25 
cluster  in  vivo,  we  analyzed  expression  of  the  miRNA  cluster  in 
transgenic  mice  in  which  Sixl  was  induced  (using  doxycycline)  in 
the  mammary  gland  (Sixl  +  Dox)  and  in  control  animals  (Ctrl  + 
Dox)^  (See  Supplementary  Figure  2  for  expression  levels  of  Sixl  in 
the  transgenic  mammary  glands)  and  found  that  all  three  miRNAs, 
miR-106b,  miR-93  and  miR-25,  are  overexpressed  in  Sixl 
transgenic  mammary  glands  as  compared  with  control  mammary 


glands  (Figure  2c).  Importantly,  Sixl  transgenic  mice  develop 
aggressive  mammary  carcinomas  that  display  multiple  histological 
subtypes  as  well  as  an  induction  of  an  EMT.^ 

The  miR-106b-25  Cluster  targets  Smad7  for  repression 
It  was  previously  shown  that  the  miR-106b-25  cluster  has  the 
ability  to  overcome  TGF-p-mediated  growth  suppression  via 
repression  of  the  cell  cycle  inhibitor  p21,  and  the  pro-apoptotic 
factor  Bim.^^  In  addition  to  the  known  role  of  the  cluster  in  TGF-p 
growth  inhibition,  using  target  prediction  analysis,  we  found  that 
the  miR-106b-25  cluster  might  also  have  a  role  in  activating  the 
TGF-p  pathway,  providing  an  attractive  mechanism  by  which  Sixl 
could  mediate  the  switch  in  TGF-p  signaling  from  tumor 
suppressive  to  tumor  promoting.  Indeed,  based  on  seed  sequence 
alignment,  the  inhibitory  Smad7  (I-Smad7)  mRNA  is  a  target  for  all 
three  miRNAs  in  the  cluster.  Smad7  antagonizes  TGF-p  signaling 
through  multiple  mechanisms,  including  binding  to  TGF-p  type  I 
receptor  (TpRI)  and  interfering  with  recruitment  and  downstream 
phosphorylation  and  activation  of  the  receptor-Smads,  Smad2  and 
Smacl3.^^  Additionally,  Smad7  also  functions  to  recruit  E3  ubiquitin 
ligases  to  TpRI,  resulting  in  its  degradation.^®  Therefore,  repression 
of  Smad7  by  the  miR-106b-25  miRNAs  would  be  expected  to 
activate  the  TGF-p-signaling  pathway,  which  is  known  to  occur 
downstream  of  Sixl.^ 

To  determine  if  Smad7  is  downregulated  in  response  to  Sixl, 
we  first  performed  qRT-PCR  on  clonal  isolates  of  MCF7-Ctrl 
and  MCF7-Six1  cells  and  demonstrated  that  Smad7  expression 
is  indeed  reduced  in  MCF7-Six1  cells,  where  the  miR-106b-25 
cluster  is  overexpressed  (Figures  3a  and  d  shows  protein  level 
differences).  To  further  determine  whether  the  cluster  of  miRNAs 
can  directly  affect  Smad7  levels,  we  generated  MCF7  cell  lines 
stably  overexpressing  the  genomic  region  of  the  cluster  (MCF7- 
Cluster),  or  control  MCF7  cells  expressing  either  the  empty  vector 
(MCF7-EV)  or  a  non-silencing  (scrambled  control)  vector  (MCF7- 
NS).  Importantly,  stable  populations  expressing  the  cluster  were 
chosen  to  overexpress  each  miRNA  in  the  cluster  only  2  to  3-fold, 
similar  to  what  is  observed  with  Sixl  overexpression  (Supplemen¬ 
tary  Figure  3).  Transfection  of  a  Smad7-3'UTR-luciferase  construct 
into  these  cell  lines  demonstrates  that  the  miR-106b-25  cluster 
inhibits  the  3'UTR  of  Smad7  (Figure  3b).  Additionally,  a  decrease  in 
Smad7  protein  in  MCF7-Cluster  cells  is  observed  when  compared 
with  I\/1CF7-EV  and  MCF7-NS  cells,  demonstrating  that  a  2-3  fold 
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Figure  1.  A  miRNA  microarray  identifies  the  miR-106b-25  cluster  family  members  as  upregulated  by  Sixl.  (a)  miRNAs  that  are  significantly  up- 
or  downregulated  (Rvalue  <0.05)  in  I\/1CF7-Six1  vs  MCF7-Ctrl  cells  as  determined  by  a  miRNA  profiling  array,  (b)  Schematic  representation  of 
the  miR-106b-25  cluster  of  miRNA  (miR-106b,  miR-93  and  miR-25)  within  the  13th  intron  of  the  MCM7  gene. 
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Figure  2.  Sixl  regulates  the  nniR-106b-25  Cluster,  (a)  Stable  overexpresson  of  Sixl  in  MCF?  cells  leads  to  an  increase  in  miR-106b,  miR-93 
and  miR-25  as  determined  using  qRT-PCR.  Data  are  represented  as  the  mean  ±s.e.m.  of  three  individual  I\/1CF7-Six1  and  I\/1CF7-Ctrl  clones, 
(b)  Knockdown  of  Sixl  in  21 PT  cells  using  Sixl -specific  siRNA  (siSixl,  50  nm  and  lOOnm)  leads  to  a  decrease  in  expression  of  all  three  miRNA  in 
the  miR-106b-25  Cluster  when  compared  with  a  control  knockdown  (siNeg).  For  qRT-PCR  analysis,  the  average  of  three  replicates  ±s.d.  is 
shown,  (c)  RNA  was  isolated  from  the  mammary  glands  of  bitransgenic  mice  in  which  Sixl  was  induced  with  doxycycline  (Sixl  +  Dox)  versus 
single  transgenic  MTB  control  mice  also  treated  with  Dox  (Ctrl  +  Dox),  but  unable  to  express  Sixl .  qRT-PCR  performed  on  the  isolated  RNA  for 
the  miR-106b-25  miRNAs  demonstrates  an  increase  in  expression  of  all  three  miRNAs  in  the  Sixl  +  Dox  mammary  glands,  which  express  high 
levels  of  the  Sixl  transgene  (Supplementary  Figure  2)  as  compared  with  Ctrl  +  Dox  control  mammary  glands.  n  =  3  mice  for  each  condition, 
and  each  miRNA  was  normalized  to  U6  RNA.  P  values  represent  statistical  analysis  using  a  paired  f-test. 
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Figure  3.  The  miR-1 06b-25  miRNAs  repress  Smad7.  (a)  qRT -  PCR  reveals 
a  decrease  in  Smad7  mRNA  in  I\/1CF7-Six1  cells  versus  I\/1CF7-Ctrl  cells. 
Data  shown  are  the  result  of  three  replicate  qRT-PCR  reactions  (±  s.d.), 
normalized  to  cyclophilin  B  mRNA  levels,  (b)  A  renilla  luciferase  reporter 
containing  the  3'UTR  of  Smad7  was  transfected  into  I\/1CF7  cells 
containing  an  EV  (I\/1CF7-EV),  a  NS  control  (I\/1CF7-NS),  or  the  miR-1 06b- 
25  cluster  (I\/1CF7-Cluster).  Measurement  of  renilla  luciferase  normalized 
to  firefly  luciferase  (present  on  the  same  vector  but  expressed  from  a 
different  promoter)  demonstrates  a  significant  repression  of  the  Smad7 
3'UTR  in  response  to  miR-1 06b-25  expression,  (c)  Western  blot  analysis 
demonstrates  that  MCF7-Cluster  cells  have  decreased  Smad7  protein  as 
compared  with  MCF7-EV  and  MCF7-NS  cells,  (d)  MCF7-Ctrl  and  MCF7- 
Sixl  cells  treated  with  miRNA  inhibitors  towards  all  three  miRNA  (ALL), 
miR-1 06b  and  miR-93,  show  a  de-repression  of  Smad7  protein  in  MCF7- 
Sixl  cells.  P  values  represent  statistical  analysis  using  a  paired  f-test. 


increase  in  the  miR-1 06b-25  cluster  can  downregulate  endogen¬ 
ous  Smad7  (Figure  3c).  Conversely,  treatment  of  MCF7-Six1  cells 
with  transient  inhibitors  against  the  individual  miRNAs  leads  to  a 
de-repression  of  Smad7  protein,  with  miR-1 06b  and  miR-93  being 
the  major  mediators  of  this  effect  (Figure  3d).  Efficacy  of  the 
miRNA  inhibitors  is  demonstrated  by  relative  activity  of  luciferase 
reporters,  containing  target  sites  for  each  miRNA,  in  the  inhibitor 
transfected  cells  (Supplementary  Figure  4). 

The  miR-1 06b-25  Cluster  activates  TGF-p  signaling 
Because  the  miR-1 06b-25  cluster  represses  Smad7  in  MCF7  cells, 
and  because  Sixl  overexpression,  which  leads  to  increased  levels 
of  the  miR-1 06b-25  miRNAs,  activates  TGF-p  signaling,^  we  asked 
whether  this  cluster  of  miRNAs,  which  is  known  to  overcome  TGF- 
p-mediated  growth  inhibition,  is  also  sufficient  to  activate  the  TGF- 
P  pathway.  Indeed,  with  both  transient  and  stable  overexpression 
of  the  miR-1 06b-25  cluster,  we  observed  an  increase  in  TpRI 
protein  levels  (Figures  4a  and  b),  as  well  as  an  increase  in  activated 
TGF-p  signaling  as  measured  by  p-Smad3  levels  (Figure  4b).  To 
further  determine  if  the  miR-1 06b-25  miRNAs  are  necessary  for  the 
previously  observed  induction  of  TGF-p  signaling  by  Sixl,  we 
utilized  a  stable  lentiviral  miRNA  knockdown  system  (miRZip)  to 
inhibit  the  miRNAs  within  the  cluster  either  individually  or 
together  (as  a  control,  we  used  a  scrambled  sequence,  miRZip- 
SCR).  Efficacy  of  the  miRZips  was  demonstrated  by  examining 
their  effects  on  endogenous  targets  of  the  miR-1 06b-25  cluster, 
p21  and  BIM  (Supplementary  Figure  5).  Inhibition  of  miR-93,  as 
well  as  the  entire  cluster  in  MCF7-Six1  cells  reverses  the  Sixl- 
induced  increase  in  TpRI  (Figure  4c),  and  inhibition  of  miR-1 06b, 
miR-93,  as  well  as  the  entire  cluster  reverses  the  Sixl -induced 
increase  in  p-Smad3  (Figure  4d). 

Interestingly,  previous  reports  have  identified  the  miR-1 06b-25 
cluster  miRNAs  as  targeting  the  TGF-p  type  II  receptor  (TpRII), 
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Figure  4.  The  nniR-106b-25  cluster  activates  TGPp  signaling,  (a)  Transient  and  (b)  stable  overexpression  of  the  miR-106b-25  cluster  in  IVICF? 
cells  leads  to  increased  expression  of  TpRI  protein  over  controls  and  an  increase  in  phosphorylated  Smad3  (p-Smad3).  (c)  I\/1CF7-Ctrl  and 
I\/1CF7-Six1  cells  expressing  stable  miRZip  inhibitors  targeting  the  miR-106b-25  miRNAs  individually  and  together  (miRZip-Cluster)  show  a 
reversal  of  the  Sixl -induced  increase  in  TpRI  protein  in  miRzip-93  and  miRZip-Cluster  treated  I\/1CF7-Six1  cells  as  compared  with  scramble 
miRZip  controls  (miRZip-SCR).  (d)  Introduction  of  miRZip-106b,  miRZip-93  and  miRZip-Cluster  into  I\/1CF7-Six1  cells  reverses  the  Sixl -induced 
increase  in  p-Smad3  levels,  without  affecting  the  Sixl -induced  increase  in  total  Smad3  levels,  (e)  A  real-time  PCR  array  containing  TGF-p 
transcriptional  targets  shows  enrichment  for  TGF-p  target  gene  expression  in  MCF 7-Cluster  cells  over  I\/1CF7-NS  cells.  Data  are  represented  as 
fold  change  expression  of  MCF  7-Cluster  compared  with  control  MCF7-NS  from  three  replicate  plates  of  each  condition.  Genes  shown  had  at 
least  a  twofold  induction  in  MCF7-Cluster  cells  as  compared  with  MCF7-NS  cells. 


resulting  in  repression  of  this  protein.^^  Analysis  of  TpRII  protein  in 
our  MCF7-Cluster  cells  did  not  show  a  repression  of  TpRII  protein 
as  compared  with  MCF7-NS  cells  (Supplementary  Figure  6). 
Similarly,  we  also  did  not  observe  significant  downregulation  of 
TpRII  protein  levels  in  MCF7-Six1  versus  MCF7-Ctrl  cells  (Supple¬ 
mentary  Figure  6).  Thus,  2-3  fold  overexpression  of  the  miR-106b- 
25  cluster  in  MCF7  cells  leads  primarily  to  alterations  in  the  TGF-p 
pathway  that  would  be  expected  to  be  activating,  as  opposed  to 
inactivating. 

To  analyze  global  changes  in  TGF-p  signaling,  we  performed 
microarray  analysis  on  MCF7-NS  versus  MCF7-Cluster  cells  to 
examine  whether  the  presence  of  the  miR-106b-25  cluster  alters 
the  TGF-p  response  signature  (TpRS),^^  similar  to  what  is  observed 
with  Sixl  overexpression.^  Hierarchical  clustering  confirmed 
differential  regulation  of  many  of  the  genes  in  the  TGF-p  response 
signature  between  MCF7-NS  and  MCF7-Cluster  cells,  demonstrat¬ 
ing  that  TGF-p  signaling  is  indeed  altered  in  response  to  miR- 
106b-25  overexpression  (Supplementary  Figure  7).  In  addition. 


we  performed  a  qRT-PCR  array  to  examine  alterations  in 
expression  of  TGF-p  target  genes  using  MCF7-NS  and  MCF7- 
Cluster  cells,  demonstrating  that  TGF-p  signaling  is  clearly 
activated  downstream  of  the  miR-106b-25  cluster,  as  numerous 
TGF-p  transcriptional  targets  are  upregulated  in  MCF7  cells 
overexpressing  the  cluster  (Figure  4e).  Of  the  84  genes  responsive 
to  TGF-p  signaling  on  the  array,  47  were  upregulated  1.5  fold  or 
more  in  MCF7-Cluster  cells,  8  genes  were  downregulated,  and  the 
rest  remained  unchanged  (Supplementary  Figure  8).  Together, 
these  data  demonstrate  that  the  miR-106b-25  cluster  is  capable  of 
activating  TGF-p  signaling  in  breast  cancer  cells,  and  suggest  that 
this  one  cluster,  which  can  also  overcome  the  growth  suppressive 
effects  of  TGF-p,  may  be  responsible  for  the  switch  in  TGF-p 
signaling  from  tumor  suppressive  to  tumor  promoting. 

The  miR-106b-25  cluster  induces  EMT-like  changes 
We  previously  reported  that  Sixl  overexpression  leads  to  an 
induction  of  EMT,  which  is  dependent  on  TGF-p  signaling.^ 
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As  the  nniR-106b-25  cluster  is  sufficient  to  activate  TGF-p  signaling, 
we  asked  if  this  cluster  is  sufficient  to  mediate  phenotypes 
associated  with  EMT.  One  of  the  hallmarks  of  EMT  is  the  loss  of 
membranous  E-cadherin  from  the  adherens  junctions.  We  thus 
analyzed  the  subcellular  localization  of  E-cadherin  from  each  cell 
line  and  demonstrated  that  E-cadherin  is  indeed  decreased  in  the 
insoluble,  or  membrane  bound,  fraction  of  the  cell  in  response  to 
miR-106b-25  overexpression  (Figure  5a).  p-Catenin,  which  is 
normally  in  the  membranous  adherens  junctions  with  E-cadherin, 
is  also  decreased  in  the  insoluble  fraction  of  MCF7-Cluster  cells 
(Figure  5a).  As  redistribution  of  p-catenin  away  from  the 
membrane  may  result  in  its  increased  nuclear  localization  and 
subsequent  ability  to  activate  transcription,  we  measured 
P-catenin  transcriptional  activity  using  the  TOP-flash  luciferase 


reporter.  Concomitant  with  the  loss  of  p-catenin  from  the 
membrane,  MCF7-Cluster  cells  also  exhibit  an  increase  in  TOP- 
flash  reporter  activity  over  I\/1CF7-NS  cells,  similar  to  the  phenotype 
observed  with  Sixl  overexpression  (Figure  5b).  To  determine  if  the 
miR-106b-25  miRNAs  are  necessary  for  the  Sixl -induced  increases 
in  p-catenin  transcriptional  activation,  we  treated  MCF7-Six1  cells 
with  inhibitors  toward  all  three  miRNA  (miR-106b,  miR-93  and 
miR-25)  together  and  measured  TOP-flash  activity.  A  repression  of 
TOP-flash  activity  in  this  context  demonstrates  that  Sixl  is 
dependent  on  the  miR-106b-25  miRNAs  to  induce  p-catenin 
transcriptional  activity  (Figure  5c).  Furthermore,  the  relocalization 
of  E-cadherin  and  p-catenin  away  from  the  membrane  in  response 
to  miR-106b-25  cluster  expression  was  confirmed  using  immuno¬ 
fluorescence  (Figure  5d). 


a 


MCF7-NS  MCF7-Cluster 


b  p-Catenin  Transcriptional  Activity  C 


p-Catenin  Transcriptional  Activity 


-  - - 

E-Cadherin 

Z) 

30 

_  _ _ 

p-Catenin 

< 

P-actin 

;> 

20 

< 

w 

10 

-  - - = 

— 

E-Cadherin 

K 

p-Catenin 

Ql 

P 

0 

MCF7-Ctrl  MCF7-Six1 


B-Catenin  Merge 


e 


□  MCF7-Ctrl  ^  MCF7-Six1  ■  MCF7-EV 
BSa  MCF7-NS  ^  MCF7-Cluster 


Figure  5.  The  miR-106b-25  cluster  mediates  features  of  EMT  (a)  miR-106b-25  overexpression  results  in  loss  of  E-cadherin  and  p-catenin  from 
the  insoluble  (cytoskeleton-associated)  protein  fraction  of  the  cell  as  determined  by  western  blot  analysis,  (b)  MCF7-Cluster  cells  show 
increased  activity  of  the  p-catenin  responsive-luciferase  reporter  TOP-Flash,  normalized  to  Renilla  luciferase  activity.  Data  are  from  three 
replicates  ±s.d.  (c)  The  Sixl -induced  increase  in  TOP-Flash  activity  is  reversed  with  inhibition  of  the  miR-106b-25  cluster  using  transient 
hairpin  inhibitors  (Dharmacon,  miRidian,  Lafayette,  CO,  USA),  (d)  Immunofluorescence  to  detect  E-cadherin  and  p-catenin  in  MCF7-NS  vs 
MCF7-Cluster  cells.  After  incubation  with  the  E-cadherin  (left  panel)  or  p-catenin  (right  panel)  antibodies,  a  FITC  conjugated  anti-secondary 
antibody  was  added  to  detect  the  signal.  Dapi  images  merged  with  FITC  are  shown  along  with  the  FITC  images  ( x  40  magnification). 

(e)  Expression  of  genes  involved  in  EMT  is  enriched  in  MCF7-Cluster  cells  versus  MCF7-NS  cells  as  shown  by  a  qPCR  array  analysis.  Gene 
expression  is  shown  by  fold  change  from  three  replicate  plates  per  group.  Genes  shown  exhibited  a  twofold  or  greater  induction  in  MCF7- 
Cluster  cells  when  compared  with  MCF7-NS  cells,  (f)  Expression  of  the  miR-106b-25  Cluster  in  MCF7  cells  results  in  decreased  adhesion  to  cell 
matrix  proteins  Collagen  I,  Collagen  IV  and  Fibronectin,  similar  to  what  is  observed  with  Sixl  overexpression.  P  values  represent  statistical 
analysis  using  a  paired  f-test  (*^0.05,  **^0.01,  ***^0.001). 
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To  deternnine  whether  additional  genes  associated  with  EMT  are 
differentially  regulated  between  MCF7-NS  and  MCF7-Cluster  cells, 
we  performed  qRT-PCR  array  analysis  on  RNA  isolated  from  each 
cell  line.  Indeed,  miR-106b-25  expression  led  to  a  2-fold  or  greater 
upregulation  of  a  number  of  genes  associated  with  EMT,  including 
Jagi,  MMP9  and  Vimentin,  amongst  others  (Figure  5e,  Supple¬ 
mentary  Figure  9  shows  all  genes  regulated  1.5-fold  or  more). 
Finally,  we  observed  a  decrease  in  cell -matrix  adhesion  to 
Collagen  I,  Collagen  IV  and  Fibronectin  in  the  MCF7-Cluster  cells 
as  compared  with  the  MCF7-NS  and  MCF7-EV  cells,  similar  to  the 
decrease  in  adhesion  observed  with  Sixl  overexpression  in  MCF7 
cells  (Figure  5f).  Together,  these  data  suggest  that  the  miR-106b- 
25  cluster  alone  can  induce  an  EMT-like  phenotype. 


The  miR-106b-25  cluster  increases  tumor-initiating  cell 
characteristics 

Many  genes  that  induce  EMT-like  phenotypes  also  induce  TIC 
phenotypes.^^  Indeed,  we  recently  demonstrated  that  Sixl 
induces  a  TIC  phenotype  in  both  a  transgenic  mouse  model^ 
and  when  overexpressed  in  MCF7  cells.^^  To  determine  if  the 
miR-106b-25  miRNAs,  which  can  induce  properties  of  EMT,  can 
also  induce  TIC  characteristics,  we  performed  flow  cytometry  for 
the  cell  surface  TIC-associated  markers  CD24  and  CD44,^^  and 
found  that  miR-106b-25  did  indeed  increase  the  percentage  of 
CD24IOW  CD44+  cells,  similar  to  Sixl  (Figure  6a).  Additionally, 
secondary  tumorsphere  assays  performed  with  MCF7-Cluster  and 
MCF7-NS  cells  demonstrated  that,  similar  to  MCF7-Six1  cells. 
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Figure  6.  The  miR-106b-25  cluster  increases  TIC  characteristics,  (a)  Overexpression  of  the  miR-106b-25  cluster  in  MCF7  cells  is  sufficient 
to  increase  the  CD24low/CD44  +  population,  similar  to  what  is  observed  with  Sixl  overexpression  in  MCF7  cells,  (b)  Expression  of  the 
miR-106b-25  cluster  in  MCF7  cells  is  sufficient  to  increase  tumorsphere  formation,  a  measurement  of  self-renewal  capability,  similar  to  what  is 
observed  with  Sixl  overexpression  in  MCF7  cells,  (c)  MCF7-Cluster  cells  transplanted  into  the  4th  mammary  fat  pad  of  NOD-SCID  mice  at 
limiting  dilutions  have  an  increased  ability  to  initiate  tumors  when  compared  with  MCF7-NS  cells,  (d)  Inhibition  of  the  miR-106b-25  cluster  in 
MCF7-Six1  cells  (MCF7-Six1 -Zip-Cluster)  reduces  the  ability  of  the  cells  to  initiate  tumors,  back  to  levels  observed  in  MCF7-Ctrl  cells  (MCF7-Ctrl- 
Zip-SCR).  The  estimated  frequency  of  tumorigenic  cells  was  calculated  by  limiting-dilution  analysis  as  described  in  the  methods,  (e)  Genes 
important  for  stem  cell  maintenance,  growth  and  differentiation  are  increased  in  MCF7-Cluster  cells  as  compared  with  MCF7-NS  cells,  as 
determined  by  a  stem  cell  qPCR  array.  Data  are  represented  as  fold  change  in  gene  expression  in  MCF7-Cluster  cells  as  compared  with  control 
MCF7-NS  cells  from  three  replicate  plates  of  each  condition.  The  Genes  shown  exhibited  at  least  a  twofold  induction  in  MCF7-Cluster  cells  as 
compared  with  the  MCF7-NS  cells. 


Oncogene  (2012),  1-10 


©  2012  Macmillan  Publishers  Limited 


The  nniR-106b-25  cluster  activates  TGF-p  Signaling 
AL  Smith  et  al 


I\/1CF7-Cluster  cells  could  increase  tumorsphere  formation,  a 
measure  of  functional  TICs  within  a  population  (Figure  6b).  To 
further  test  for  functional  TIC  characteristics  using  an  in  vivo  assay, 
we  injected  cells  at  limiting  dilutions  into  the  mammary  fat  pad  of 
NOD-SCID  mice.  I\/1CF7-Cluster  cells  were  able  to  initiate  tumors 
with  a  greater  frequency  than  MCF7-NS  cells,  both  when  1000  and 
100  cells  were  injected  (Figure  6c).  In  order  to  determine  if  the 
miR-106b-25  cluster  was  necessary  for  Sixl -induced  increases  in 
TICs  in  vivo,  we  utilized  our  MCF7-Six1-miRZip-Cluster  cells  (in 
which  all  three  miRNAs  are  inhibited),  and  transplanted  these  cells 
into  the  mammary  gland  at  limiting  dilutions,  along  with  miRZip- 
SCR  controls  in  both  MCF7-Ctrl  and  MCF7-Six1  cells.  Inhibition  of 
the  miR-106b-25  miRNAs  in  MCF7-Six1  cells  demonstrates  a 
reduction  in  tumor  initiating  ability  back  to  similar  levels  of  TIC 
frequency  seen  in  MCFZ-Ctrl  cells  (Figure  6d).  Lastly,  we  performed 
a  human  stem  cell  qRT-PCR  array.  Of  the  84  genes  on  this  array, 
53  were  upregulated  more  than  1.5  fold  in  MCF7-Cluster  cells 
(Figure  6e  and  Supplementary  Figure  10).  Together  these  data 
demonstrate  for  the  first  time  a  role  for  the  miR-106b-25  cluster  in 
both  EMT  and  in  increased  TIC  capacity. 
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The  miR-106b-25  cluster  correlates  with  Six,  activated  TGF-p 
signaling,  and  shortened  time  to  relapse  in  human  breast  cancer 

To  determine  if  the  Six1/miR-106b-25/activated  TGF-p  signaling 
axis  is  relevant  to  human  breast  cancer,  we  obtained  human 
breast  cancer  tissue  arrays  containing  71  cases  of  invasive  ductal 
carcinoma,  with  matched  cases  on  which  we  had  previously 
performed  immunohistochemistry  using  an  antibody  generated 
against  human  Sixl  and  human  Smad3  (nuclear  Smad3  was 
scored  as  an  indicator  of  activated  TGF-p  signaling).^  We  next 
performed  in  situ  hybridization  for  miR-106b  (as  a  representative 
of  the  miR-106b-25  cluster),  to  compare  expression  of  this  cluster 
family  member  with  nuclear  staining  obtained  with  the  Sixl 
antibody  and  with  activated  TGF-p  signaling.  Importantly,  miR- 
106b  and  Sixl  significantly  correlate  in  breast  cancer  tissues 
{P  =  0.0028  Spearman's  /?  =  0.3927)  (Figure  7a),  as  do  miR-106b 
and  nuclear  Smad3  (P  =  0.001 7,  Spearman's  /?  =  0.3972) 
(Figure  7b).  In  addition,  the  greatest  percentage  of  tumors 
exhibited  activated  TGF-p  signaling  when  both  miR-106b  and 
Sixl  were  highly  expressed  (64.7%  show  increased  nuclear  Smad3 
when  both  miR-106b  and  Sixl  are  high)  (Figure  7c).  Of  note. 


a 


Six1  correlates  with  miR-106b 
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Figure  7.  The  miR-106b-25  cluster  correlates  with  Sixl  expression,  activated  TGFp  signaling  and  shortened  time  to  relapse  in  human  breast 
cancers.  Human  breast  cancer  tissue  arrays  were  previously  immunostained  with  an  anti-Sixl  antibody  (Atlas)  and  a  Smad3  antibody  (Zymed) 
as  previously  described^  and  the  staining  was  scored  for  nuclear  Sixl  and  Smad3  on  a  scale  of  0-4.  A  serial  section  array  was  also  stained  for 
miR-106b  expression  by  in  situ  hybridization.  Expression  of  miR-106b  was  scored  on  a  scale  of  0-4  and  compared  with  Sixl  and  nuclear 
Smad3  scores  in  the  same  tissues,  (a)  Results  show  that  miR-106b  and  Sixl  correlate  in  human  breast  cancers  (b)  as  do  miR-106b  and  nuclear 
Smad3.  (c)  When  the  expression  of  all  three  molecules  are  considered,  the  highest  percentage  of  nuclear  Smad3  can  be  found  when  both  Sixl 
and  miR-106b  are  highly  expressed.  P-values  obtained  using  Spearman's  correlation  analysis,  (d)  In  a  miRNA  expression  dataset  of  216  early- 
invasive  breast  cancers,^^  patients  whose  tumors  express  both  high  miR-106b  and  high  miR-93  show  a  significantly  reduced  time  to  relapse. 
The  median  value  for  miR-106b  and  miR-93  was  used  to  divide  the  samples  into  high  (above  median)  and  low  (below  median)  miRNA 
expression.  P  value  was  calculated  by  log-rank  analysis.  A  full  colour  version  of  this  figure  is  available  at  the  Oncogene  journal  online. 
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the  Sixl  antibody  used  in  these  experiments  was  generated  against 
a  region  of  Sixl  that  may  allow  cross-reactivity  with  other  Six  family 
members,  as  they  are  highly  conserved.  Therefore,  we  can  only 
confidently  say  that  Six  family  expression  correlates  with  miR-106b 
in  human  tumors.  However,  examination  of  the  van  de  vijver  public 
breast  cancer  gene  expression  dataset^^  suggests  that  Sixl  is  more 
abundantly  expressed  than  other  Six  family  members  in  human 
breast  tumors  (data  not  shown).  Together,  these  data  strongly 
suggest  a  critical  role  for  miR-106b-25  in  the  Sixl -induced  activation 
of  TGF-p  signaling  in  human  breast  cancer. 

To  explore  the  prognostic  value  of  these  miRNA  in  human 
breast  cancers,  we  examined  a  publicly  available  dataset 
comprised  of  miRNA  expression  in  early-invasive  breast  cancers.^^ 
Figure  7d  demonstrates  that  patients  whose  tumors  express  high 
miR-106b  and  high  miR-93  together  have  a  significantly  shortened 
time  to  relapse.  Analysis  of  individual  miRNA  expression  in  these 
tumors  also  demonstrates  a  significant  correlation  with  high 
miR-93  (Supplementary  Figure  11B),  as  well  as  a  trend  toward 
shortened  time  to  relapse  with  high  miR-106b  and  high 
expression  of  all  three  miRNA  (Supplementary  Figure  11A,D). 
However,  miR-25  expression  does  NOT  demonstrate  any  differ¬ 
ence  in  patient  outcome  (Supplementary  figure  11C),  further 
suggesting  that  miR-106b  and  miR-93  are  the  primary  regulators 
of  this  response. 

To  further  explore  the  role  of  miR-106b  and  miR-93  in  human 
breast  cancers,  we  also  utilized  gene  set  enrichment  analysis 
(GSEA)  to  analyze  breast  tumors  that  express  high  levels  of  miR- 
106b  and  miR-93.  Of  the  enriched  gene  sets,  we  found  high  miR- 
106b  and  high  miR-93  to  significantly  correlate  with  stem  cell 
related  gene  sets^^  as  well  as  EMT  enriched  gene  sets,^^  further 
supporting  an  important  role  for  these  miRNAs  in  mediating  stem 
cell  and  EMT-like  processes  in  human  breast  cancer  (Supplemen¬ 
tary  Figure  12). 


DISCUSSION 

In  the  present  study,  we  identify  a  cluster  of  miRNA  previously 
shown  to  overcome  TGF-p-mediated  growth  inhibition,^  ^  the  miR- 
106b-25  cluster,  as  a  target  of  Sixl.  The  data  herein  not  only 
provide  a  mechanism  for  how  Sixl  may  silence  TGF-p-mediated 
growth  inhibition  in  breast  cancers,  but  also  for  how  Sixl  activates 
the  TGF-p  pathway.  This  is  the  first  demonstration  that  the  same 
miRNA  cluster  that  overcomes  TGF-p-mediated  growth  suppres¬ 
sion  can  in  fact  also  promote  TGF-p  signaling. 

Indeed,  we  show  that  miR-106b-25  miRNAs  can  target  the 
TGF-p  inhibitor  Smad7,  and  that  upregulation  of  miR-106b-25 
leads  to  an  increase  in  TpRI.  This  increase  in  TpRI  likely  occurs,  at 
least  in  part,  due  to  Smad7  downregulation,  as  Smad7  is  known  to 
mediate  degradation  of  the  TpRI  protein. Recently,  we  demon¬ 
strated  that  upregulation  of  TpRI  protein  is  necessary  and 
sufficient  for  TGF-p  activation  and  the  induction  of  EMT  down¬ 
stream  of  Sixl  in  MCF7  cells.^  Consistent  with  our  previous  data 
demonstrating  that  TpRI  overexpression  is  sufficient  for  TGF-p 
pathway  activation,  we  also  observe  an  activation  of  this  pathway 
with  overexpression  of  the  miR-106b-25  miRNAs  alone  as 
evidenced  by  an  increase  in  transcriptional  targets  of  the  TGF-p 
signaling  pathway  (Figure  4e). 

Several  lines  of  evidence  have  demonstrated  that  the  miR-106b- 
25  cluster  and  its  individual  miRNAs  have  pro-oncogenic 
functions,  including  mediating  pro-proliferative  and  anti-apoptotic 
phenotypes.^ Our  results  in  this  study  expand  the  oncogenic 
potential  of  this  miRNA  cluster  by  demonstrating  for  the  first  time 
that  these  miRNA  can  also  induce  properties  of  EMT  and  TIC 
characteristics.  The  EMT  changes  induced  by  these  miRNA  are 
consistent  with  the  oncogenic  EMT  phenotype  induced  by  Sixl  in 
MCF7  cells  (Figure  5),  suggesting  that  overexpression  of  the  miR- 
106b-25  miRNAs  may  partly  contribute  to  the  induction  of  EMT 
downstream  of  Sixl . 


It  is  well  recognized  that  the  induction  of  EMT  leads  to  an 
increase  in  stem/progenitor  cell  properties.^^  Indeed,  Sixl 
transgenic  mice  whose  mammary  tumors  display  features  of 
EMT,  also  demonstrate  an  increase  in  the  stem/progenitor  cell 
population.^  Increased  expression  of  the  miR-106b-25  miRNAs  in 
the  mammary  glands  of  Sixl  transgenic  mice  suggests  a  possible 
role  for  these  miRNA  in  regulation  of  the  stem/progenitor  pool 
(Figure  2c).  Our  results  show  for  the  first  time  that  the  miR-106b- 
25  miRNAs  are  sufficient  to  increase  TIC  capacity,  and  that  they  are 
required  for  the  ability  of  Sixl  to  induce  TIC  characteristics 
in  vivo.  Interestingly,  many  recent  studies  have  implicated  the  miR- 
106b-25  cluster  in  stem/progenitor  cell  biology.^^'^®'^^  Of  interest, 
it  was  recently  demonstrated  that  miR-106b  and  miR-93,  along 
with  iPSC  transcription  factors,  can  enhance  reprogramming  of 
somatic  cells  into  induced  pluripotent  stem  cells.^^  The  iPSC 
phenotype,  however,  is  dependent  on  downregulation  of  TGF-p 
signaling,  where  miR-106b  and  miR-93  target  TpRII  (in  a  mouse 
embryonic  fibroblast  context),  while  the  TIC  phenotype  is  known 
to  be  associated  with  an  upregulation  of  TGF-p  signaling.^^  In  this 
study,  we  did  not  observe  a  downregulation  of  TpRII  in  response 
to  Sixl  or  miR-106b-25  overexpression  in  MCF7  cells  (Supplemen¬ 
tary  Figure  6).^  However,  interestingly,  both  the  iPSC  data  and  our 
data  suggest  a  role  for  the  miR-106b-25  cluster  in  the  induction 
of  stem  cell  properties,  possibly  through  their  ability  to  regulate 
TGF-p  signaling  in  a  context  dependent,  and  seemingly  opposite, 
manner. 

As  the  increase  of  TIC  capacity  is  known  to  be  associated  with 
tumor  recurrence,^^  our  data  are  further  strengthened  by  our 
analysis  of  miRNA  microarray  data,  which  demonstrate  that  high 
miR-106b  and  high  miR-93  expression  significantly  correlate  with  a 
shortened  time  to  relapse  in  human  breast  cancer  patients 
(Figure  7d).  Of  note,  consistent  with  miRNA-induced  iPSC  data,^^ 
our  data  also  suggests  that  miR-106b  and  miR-93  (which  are  in  the 
same  family  of  miRNA)  are  the  primary  inducers  of  the  TIC 
phenotype  as  there  is  no  difference  in  tumor  recurrence  when 
only  miR-25  is  highly  expressed  (Supplementary  Figure  11). 
Similarly,  in  breast  tumors  where  both  miR-106b  and  miR-93  are 
highly  expressed  together,  an  enrichment  in  stem  cell  and  EMT 
gene  signatures  can  be  observed  (Supplementary  Figure  12), 
further  strengthening  the  argument  that  the  miR-106b  family  of 
miRNAs  may  be  important  regulators  of  EMT  and  TIC  phenotypes 
in  human  breast  cancer. 

Lastly,  our  data  also  demonstrate  a  significant  correlation 
between  miR-106b,  Sixl  and  activated  TGF-p  signaling  (nuclear 
Smad3)  in  human  breast  cancers.  Critically,  we  show  that  tumors 
that  express  both  high  Sixl  and  high  miR-106b  have  the  highest 
percentage  of  activated  TGF-p  signaling  (64.7%)  (Figure  7).  Of 
note,  our  data  also  show  that  in  the  presence  of  high  Sixl  and  low 
miR-106b  expression,  23.5%  of  tumors  have  activated  TGF-p 
signaling,  as  opposed  to  5.9%  in  tumors  with  low  levels  of  Sixl. 
These  data  suggest  again  that  Sixl  may  activate  TGF-p  signaling 
through  multiple  mechanisms,  however,  the  marked  increase  in 
TGF-p  signaling  when  Sixl  and  miR-106b  are  both  highly 
expressed  strongly  suggests  that  the  Six1/miR-106b-25  axis  is 
critical  for  activation  of  TGF-p  signaling  in  human  breast  cancer. 

In  closing,  these  data  have  important  ramifications  for  breast 
cancer  treatment.  Due  to  the  traditional  difficulties  in  targeting 
transcription  factors,  such  as  Sixl,  and  the  increasing  promise  for 
miRNAs  as  therapeutic  targets,^^  the  miR-106b-25  cluster  could 
prove  to  be  an  effective  target  in  cancers  that  express  high  levels 
of  Sixl .  Furthermore,  as  TGF-p  signaling  can  be  tumor  suppressive 
or  tumor  promotional,  depending  on  context,  one  of  the  greatest 
concerns  surrounding  the  use  of  TGF-p  inhibitors  in  cancer  is  how 
to  predict  which  patients  will  benefit  in  clinical  trials.  Our  studies 
suggest  a  mechanism  that  we  hypothesize  provides  a  novel 
molecular  explanation  for  the  TGF-p  paradox  in  breast  cancers. 
Namely,  examining  breast  tumors  for  Sixl  and/or  miR-106b-25 
expression  may  ultimately  provide  a  means  to  distinguish  patients 
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likely  to  benefit  from  TGF-p  inhibitors  from  those  who  may 
actually  be  harmed  by  such  treatments. 


MATERIALS  AND  METHODS 

microRNA  microarray  and  mRNA  microarray 

Total  RNA  preparations  using  TRIzol  were  sent  to  Thermo  Scientific 
(Waltham,  MA,  USA)  for  miRNA  microarray  analysis.  Three  clonal  isolates  of 
IV1CF7-Ctrl  and  IVICF7-Six1  RNA  samples  were  submitted  in  replicates  of 
four.  miRNAs  in  which  the  difference  between  the  groups  (IVICF7-Ctrl  and 
IVICF7-Six1)  was  statistically  significant  (P-value  <0.05)  were  chosen  for 
follow-up  analysis.  For  mRNA  microarray  analysis,  RNA  was  prepared  in  the 
same  way  as  above,  and  submitted  to  the  Genomics  and  Microarray  Core 
at  the  University  of  Colorado  Denver.  The  array  was  performed  using  the 
Affymetrix  Fluman  Gene  1.0  ST  Array  (Affymetrix,  Santa  Clara,  CA,  USA). 
Analysis  of  microarray  results  have  been  described  previously.^ 

Cell  culture  and  constructs 

Generation  of  IVICF7-Ctrl  and  I\/ICF7-Six1  cell  lines  was  described 
previously.^^  The  MCF 7-Cluster  cells  were  generated  by  inserting  the 
cloned  genomic  region  of  the  miR-106b-25  cluster  into  the  miR-Express 
vector  (Open  Biosystems,  Lafayette,  CO,  USA).  EV  and  non-silencing  (NS) 
miR-express  constructs  were  obtained  from  open  biosystems.  The  miRZip- 
cluster  construct  contains  all  three  miRNA  inhibitors  (miRZip-106b,  miRZip- 
93  and  miRZip-25)  in  the  same  vector.  I\/ICF7-Ctrl  and  IVICF7-Six1  cells  were 
infected  with  miRZip  lentivirus  (System  Biosciences,  Mountain  View,  CA, 
USA).  All  MCF7  cells  were  selected  with  2.5|ig/ml  puromycin. 

Real-time  PCR 

Total  RNA  was  extracted  with  the  miRNeasy  RNA  isolation  kit  (Qiagen, 
Germantown,  MD,  USA)  following  the  manufacturer's  protocol.  For  miRNA 
quantitative  analysis,  RNA  was  reverse  transcribed  using  the  miscript  system 
(Qiagen),  and  qPCR  was  performed  with  miscript  miRNA  primers  (Qiagen).  All 
miRNA  assays  were  done  using  ssoFast  Evagreen  supermix  (Biorad,  Hercules, 
CA,  USA).  For  mRNA  qPCR,  cDNA  synthesis  was  done  with  iscript  (BioRad). 
Sixl,  Smad7  and  PPIB  primers  were  part  of  the  gene  expression  assay 
collection  (Applied  Biosystems,  Carlsbad,  CA,  USA).  All  qPCR  was  performed 
with  the  BioRad  CFX96.  Real-time  PCR  arrays  were  acquired  from 
SABiosciences  (Frederick,  MD,  USA)  (RT^  Profiler  PCR  array,  PAHS-235,  PAHS- 
405  and  PAHS-090),  and  performed  according  to  manufacturer  instructions. 

Western  blot  and  immufluorescence 

For  western  blots,  whole  cell  lysates  were  isolated  using  RIPA  buffer  as 
previously  described.^^  Antibodies  used  include:  E-cadherin  (BD  Bios¬ 
ciences,  Franklin  Lakes,  NJ,  USA),  p-catenin  (BD  Biosciences),  T^RI  (SCBT, 
Santa  Cruz,  CA,  USA),  T^RII  (Cell  Signaling,  Danvers,  MA,  USA),  p-Smad3 
(Cell  Signaling),  total  Smad3  (Invitrogen,  Grand  Island,  NY,  USA),  Bim  (Cell 
Signaling),  p21  (Cell  Signaling),  p-Actin  (sigma-aldrich,  St  Louis,  MO,  USA) 
and  p-Tubulin  (Invitrogen).  Cell  fractionation  was  performed  as  previously 
described.^^  Immunofluorescence  was  performed  as  previously  described^ 
with  E-cadherin  and  p-catenin  antibodies  above. 

Cell  adhesion  assay 

Cells  were  plated  in  96-well  plates  coated  with  Collagen  I,  Collagen  IV, 
Laminin  or  Fibronectin  (BD  Biocoat,  BD  Biosciences),  and  assays  were 
carried  out  as  previously  described.^ 
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degrees  below  the  RNA  Tm.  Detection  was  achieved  with  BM  purple  (Roche, 
Basel,  Switzerland)  solution,  and  slides  were  counterstained  in  Nuclear  Fast 
Red  (Poly  Scientific,  Bay  Shore,  NY,  USA).  Slides  were  scored  on  a  0-4  scale 
with  4  representing  the  most  intense  staining.  Scores  were  assigned 
independently  by  three  individuals  in  a  blinded  manner  and  averaged.  Serial 
sections  of  tumor  arrays  were  previously  stained  and  scored  by  a  pathologist 
on  a  scale  of  0-4  for  nuclear  Sixl  (1:100;  Atlas  antibodies,  Stockholm, 
Sweden)  and  nuclear  Smad3  (5  pg/ml;  Zymed,  San  Francisco,  CA,  USA),  using 
immunohistochemistry  protocols  previously  described.^ 


Luciferase  assays 

The  3'UTR  of  Smad7  was  cloned  into  the  psi-Check2  luciferase  reporter 
(Promega,  Madison,  Wl,  USA).  For  perfect  target  luciferase  assays,  the  exact 
complement  sequence  of  each  miRNA  was  also  cloned  into  the  psi-Check2 
luciferase  reporter.  Constructs  were  transfected  using  Lipofectamine  2000 
transfection  reagent  (Invitrogen),  and  at  48  h  lysates  were  prepared  and 
analyzed  using  the  dual  luciferase  assay  following  the  manufacturer's 
protocol  (Promega).  TOP-flash  reporter  assays  were  done  as  described 
previously^  All  luciferase  assays  were  analyzed  on  the  Modulus  Microplate 
reader  (Turner  Biosystems,  Madison,  Wl,  USA). 

TIC  assays 

Flow  cytometry  analysis  and  tumorshpere  formation  assays  were  performed 
as  described  previously^^  For  in  vivo  tumor  initiation  assays,  cells  were 
counted  and  serially  diluted  in  100  pi  of  1:1  PBS/Matrigel  (no.  354234,  BD 
Biosciences).  Diluted  cells  were  injected  underneath  the  nipple  of  the  number 
4  mammary  fat  pad  of  6-week-old  female  NOD/SCID  mice.  Tumor  formation 
was  monitored  weekly  by  palpation.  The  frequency  of  tumorigenic  cells  was 
estimated  at  a  setting  of  95%  confidence  using  the  ELDA  software  application 
at  http://bioinf.wehi.edu.au/software/elda/.  All  animal  studies  were  performed 
according  to  protocols  reviewed  and  approved  by  the  Institutional  Animal 
Care  and  Use  Committee  at  the  University  of  Colorado  Denver. 


Analysis  of  nniRNA  microarray  data 

miRNA  expression  and  clinical  outcome  data  were  acquired  from  a 
publically  available  dataset  of  216  early-invasive  primary  breast  cancers^^ 
(GSE22220).  All  samples  were  median-centered  for  each  miRNA  expression, 
and  denoted  high  expression  if  above  the  median  and  low  expression  if 
below  the  median.  Kaplan -Meier  curves  for  tumors  expressing  both  high 
miR-106b  and  high  miR-93  were  generated  using  WinStat  for  Excel  (R  Fitch 
Software,  A-Prompt  Corp,  Whitehall,  PA,  USA).  Gene  set  enrichment 
analysis  methods  are  included  in  the  Supplementary  Information. 
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In-situ  hybridization  and  Immunohistochemistry 
Breast  cancer  tissue  array  (BRC71 1)  was  purchased  from  USA,  Biomax  Inc. 
(Rockville,  MD,  USA).  In  situ  hybridization  was  performed  with  double  DIG 
labeled  miRNA  LNA  probes  (Exiqon,  Vedbaek,  Denmark)  using  manufac¬ 
turer's  one-day  protocol.^^  Modifications  include  15  pg/ml  Proteinase  K  for 
8  min,  overnight  hybridization  with  40  nM  of  LNA  probe,  and  a  formamide 
containing  hybridization  buffer  (50%  Formamide,  5x  SSC,  0.1%  Tween, 
50  pg/ml  Heparin,  500  pg/ml  yeast  tRNA)  at  hybridization  temperatures  30 
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